WM-97-105 

CFNUL-97-01 

JLAB-THY-97-09 

Covariant equations for the three-body bound state 



Alfred Stadler 1 ' 2 , Franz Gross 1 ' 3 , and Michael Frank 4 
1 College of William and Mary, Williamsburg, Virginia 23185 
2 Centro de Fisica Nuclear da Universidade de Lisboa, 1699 Lisboa Codex, Portugal? 
3 Thomas Jefferson National Accelerator Facility, Newport News, VA 23606 
^ Institute for Nuclear Theory, University of Washington, Seattle WA 98195 



Abstract 



The covariant spectator (or Gross) equations for the bound state of three 
identical spin 1/2 particles, in which two of the three interacting particles are 
always on shell, are developed and reduced to a form suitable for numerical 
solution. The equations are first written in operator form and compared 
to the Bethe-Salpeter equation, then expanded into plane wave momentum 
states, and finally expanded into partial waves using the three-body helicity 
formalism first introduced by Wick. In order to solve the equations, the two- 
body scattering amplitudes must be boosted from the overall three-body rest 
frame to their individual two-body rest frames, and all effects which arise from 
these boosts, including the Wigner rotations and p-spin decomposition of the 
off-shell particle, are treated exactly. In their final form, the equations reduce 
to a coupled set of Faddeev-like double integral equations with additional 
channels arising from the negative p-spin states of the off-shell particle. 

21.45+v, 11. 10. St, 21.10.Dr 
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I. INTRODUCTION AND OVERVIEW 



The three-body spectator (or Gross) equations were first introduced and applied to scalar 
particles in 1982 This original paper included a treatment of non-identical particles and 
an introductory discussion of the definition and role of three-body forces in a relativistic 
context. Shortly afterward, in lectures given at the University of Hannover ||, the equations 
for three identical spin 1/2 particles were written down, but many details needed for a 
practical solution of the equations were never worked out. In this paper we complete the 
development by expanding the amplitudes into partial waves and reducing the equations to 
a compact form suitable for numerical solution. The development is carried out only for 
the case when the three-body scattering amplitude can be obtained by iterating successive 
two-body interactions, so that the three-body forces of relativistic origin discussed in the 
original paper |J are neglected. However, because our covariant equations include the 
negative energy part of the Dirac propagator of the off-shell nucleon, many contributions 
are automatically included which would arise from three-body forces in a nonrelativistic 
context. 

The bound state equations we present in this paper have already been solved numeri- 
cally for a variety of cases, and some results have already been published |j|,f||. From this 
experience we know that the general development presented here is a suitable basis for a 
practical solution of the covariant three-body problem. 

In the remainder of this section we give a brief summary of the current status of nonrel- 
ativistic calculations of the binding energy of the three nucleon bound state, and a review of 
previous work on the relativistic three-body problem. Then we will give a brief summary of 
the physics underlying our spectator equations, and present the final equations. The deriva- 
tion of these results is found in the subsequent sections. In Sec. II we begin the development 
by writing the three-body equations in an operator form which is independent of the basis 
states used to describe the three-body system. In Sec. Ill we introduce basis states and 
write the equations in momentum space. In this representation the physical content of the 
equations is clear, but the equations are not in a form most convenient for numerical solu- 
tion. To solve the equations numerically it is convenient to use a partial wave decomposition 
based on the helicity states originally introduced (in a three-body context) by Wick |J and 
this is developed in detail in Sec. IV. The evaluation of the permutation operator, which 
interchanges particles between interactions and permits us to express the equation in terms 
of only one amplitude, is discussed in detail in Sec. V, and all of the results are collected 
together and the final equations given in Sec. VI. There are three appendices which discuss 
some points in detail. 



A. Brief history of the three-body bound state problem 

The first realistic nonrelativistic calculations of the triton binding energy were completed 
in the 1970's ||. Later it was shown that different methods arrived at the same results, 
and that the binding energy could be calculated to a numerical accuracy of a few keV by 
considering all nucleon- nucleon (NN) partial waves up to j = 4 fjj. Today, if three-body 
forces (3BFs) are not considered, a discrepancy of about 0.5-1.0 MeV remains between the 



3 



experimentally observed value of —8.48 MeV and values obtained from realistic nonrelativis- 
tic NN potentials. Calculations of the contribution of the A resonance to the 3BF find that 
the net effect of the A is small 0,0]. State-of-the art calculations often include in addition 
also 3BFs based on meson-nucleon interaction processes other than A excitation ||10|| . When 
the strength of phenomeno logical 3BFs is adjusted to give the correct triton binding energy, 
an excellent value is also obtained for the 4 He binding energy (and to a lesser extent other 
light nuclei up to A ~ 7) |QJ . 

However, relativistic effects should make a contribution to the binding energy at the 
level of several hundred keV. Using a mean momentum of about 200 MeV (consistent with 
nonrelativistic estimates) we expect to see corrections of the order of (v/c) 2 ~ (p/m) 2 ~ 4%. 
If this is 4% of the binding energy, then it amounts to about 300 keV. However, if relativity 
has a greater effect on the attractive a exchange part of the force (as it does in nuclear 
matter calculations using the Walecka model H|) then we might obtain an effect 10 times 
larger. 

Interest in relativistic three-body equations goes back to 1965, when Alessandrini and 
Omnes Jl3| used the Blankenbecler-Sugar equation |L4| to describe the scattering of three 
particles, and Basdevant and Kreps |H| applied their ideas to a description of the three 

II| discussed the application of the Bethe-Salpeter equation |]TTJ to 
1966. In 1968 Aaron, Amado, and Young |l£j introduced three- 
body scattering equations in which all the particles were on shell. Later, Garcilazo and his 
collaborators |19[ treated three-body bound states using the Blankenbecler-Sugar equation, 
and Garcilazo |20| applied Wick's helicity formalism to the three-body problem, and used 
it to treat the ttNN system relativistically . Recently, the size of relativistic effects were 
estimated by Rupp and Tjon [^2] using a separable kernel in the Bethe-Salpeter equation, by 
Sammarruca, Xu, and Machleidt [2~3"| using minimal relativity and the Blankenbecler-Sugar 
equation, and by the Urbana group [24J using the Schrodinger equation with corrections 
of first order in (v/c) 2 . All of these calculations include some contributions coming from 
relativistic kinematics, but none treats the full Dirac structure of the nucleons, or investigates 
effects which might arise from a realistic relativistic treatment of the NN dynamics. 



pion system. Taylor 
three-body systems in 



B. The physics behind the spectator equations 

In the absence of three-body forces, the three-body scattering amplitude (and the three- 
body bound state vertex functions) can be obtained by summing all successive two-body 
scatterings, as shown diagrammatically in Fig. 1. This summation can be organized into 
Faddeev-like equations, shown diagrammatically in Fig. 2. When the three particles are iden- 
tical, the different Faddeev subamplitudes can be obtained from each other by interchange 
of variables, leading to a single equation for a single subamplitude represented diagram- 
matically in Fig. 3. It is necessary to know the two-body scattering amplitude before the 
equation shown in Fig. 3 can be solved. More specifically, the two-body scattering amplitude 
must be known in the rest frame of the three-body system (or any other frame independent 
of the internal variables). 

The two-body amplitude is usually calculated in its own rest frame, so it must be boosted 
to the three-body rest frame before it can be used in the Faddeev equations. The velocity 
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FIG. 1. Diagrams from the infinite class of successive two-body scatterings (represented by 
the ovals) which contribute to scattering of the three-body system. We use the convention that 
the initial state is on the right and the final state on the left in each diagram. The x labels 
internal spectators, which are put on-shell in the spectator formalism. In this example, all of the 
diagrams but the first contribute to the subamplitude T 13 where particle 1 (on the top) is the last 
spectator and particle 3 (on the bottom) is the first spectator. The first diagram contributes to 
the subamplitude 7" . 

of this boost depends on the momentum of the spectator, which is one of the dynamical 
variables of the problem, and hence the boost must be known for all velocities. In the non- 
relativistic case this is trivial because the two-body amplitude is invariant under Galilean 
boosts. However, in the relativistic case this may present a problem, depending on the type 
of formalism used. Here, for the purposes of discussion, we distinguish two fundamentally 
different ways to approach relativistic calculations. In one approach, which will be referred 
to as Hamiltonian Dynamics (including light-cone methods) [^5| , some of the Poincare gen- 
erators include the interactions, and either boosts or rotations cannot be carried out exactly. 
In this method one must treat relativistic effects approximately. In a second method, which 
we will refer to as Manifestly Covariant Dynamics |[26|| , the generators are all kinematic, 
and the boosts can be done exactly. The spectator equations developed in this paper are an 
example of the latter method; we will reduce the three-body equations to a practical form 
by exploiting our ability to boost the two-body amplitudes to their rest frame. 

Of the methods discussed in the previous subsection, only the Bethe-Salpeter (BS) for- 
malism shares the property that the two-body amplitudes can be boosted exactly to their 
rest frame, and we will therefore compare the spectator equations with the corresponding 
BS equations. Both approaches conserve total four-momentum. This leaves an integra- 
tion over all independent internal four-momenta, which are two for the three-body problem. 
The three particle Bethe-Salpeter equation does not restrict any of these eight independent 
components, and after a partial wave decomposition there still remain four integrations, lead- 
ing to coupled /owr-dimensional Faddeev equations. Furthermore, these equations contain 
singularities arising from the indefinite nature of the Minkowsky metric. In the specta- 
tor formalism the two time components of the internal four-momenta are eliminated (or, 
more precisely, expressed in terms of the other variables) by requiring that two of the three 
particles be fixed to their positive energy mass shell. This reduces the number of inde- 
pendent variables to only six, and after a partial wave decomposition one obtains coupled 
iu>o-dimensional equations with a Faddeev structure. The three-body spectator equations 
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FIG. 2. Diagrammatic representation of the Faddeev equations for the amplitudes T 1% . Note 
that the spectator is identified by the solid dot. 

therefore have the same structure as nonrelativistic equations, and this is one of their most 
significant advantages. 

A particle is put on-shell when it is a spectator to the interaction of two other particles. 
When this is done systematically, two of the three particles are always on-shell. The particle 
which is off shell is the (unique) particle which has just interacted and is about to interact 
again (in a topological, not time-ordered, sense), as illustrated in Fig. 1. 

It is natural to assume that restricting particles to their mass-shell represents an approx- 
imation to the BS equation, but it can be shown that it is equivalent to a reorganization of 
the perturbation series of all ladder and crossed ladder diagrams which, in some cases, sums 
these diagrams more efficiently P|J27[]. 

In summary, the spectator equation is used because: 

(i) it sums the infinite series of all ladder and crossed ladder interactions efficiently, 

(ii) it reduces the number of independent variables to a minimum, making the covariant 




FIG. 3. Diagrammatic representation of the bound state spectator equation for three identical 
particles. Spectators are identified by the solid dots, and on-shell particles by the x. Note the 
interchange of particles 1 and 2. 
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three-body problem tractable, and 



(iii) it permits us to boost the two-body amplitudes to their rest frame and calculate 
relativistic effects exactly. 

Before we turn to the details of the derivation of the spectator equations, we present the 
equations in their final form in the next subsection. 



In the absence of 3BFs the three-body scattering amplitude is obtained from a sum 
of all successive two-body scatterings. Because the three particles are identical, each two- 
body scattering differs from the others only by a permutation, and they can therefore all be 
summed by one operator equation of the form 



where IT 1 ) is a vertex function describing the contribution to the bound state from all 
processes in which the 23 pair was the last to interact (with particle 1 a spectator), the 
two-body amplitude M 1 describes the scattering of the 23 pair, G 1 is the propagator for 
the 23 pair, and Vu is a permutation operator interchanging particles 1 and 2. These are 
labeled in Fig. 3. The factor of 2 comes from the contribution of V13 which equals the one 
of Vi2- The permutation operator rearranges the particles so that the same equation sums 
up the scattering of all pairs: 12, 23, and 13. 

The three-body spectator equations have the same structure as but incorporate 

the additional feature that the spectator is restricted to its positive energy mass-shell in 
all intermediate states. With the conventions implied above, consistency also requires that 
particle 2 be on-shell, so that two particles are always on-shell. As already stated above, we 
think of these constraints as a reorganization of Eq. (|1.1|) which will, in some cases, improve 
its convergence. The constraints are manifestly covariant, and lead to the following equation 



where the lower index labels the second on-shell particle. Hence only particle 3, the (unique) 
particle which has just left one interaction and is about to enter another one, is off-shell in 



using three-particle states. Both p-spin states (where p = + is the u spinor positive energy 
state and p = — is the v spinor negative energy state) of the off-shell particle must be 
treated. First we reduce the equation using states with definite particle helicities, similar 
to those defined by Wick ||. These three-body states will be written in the abbreviated 
form I Jl(23)p), where J is the total angular momentum of the state, p the p-spin of the 
off-shell particle, 1 = {q,Xi} (where q and Ai are the magnitude of the three-momentum 
and the helicity of the spectator in the three-body cm.), and (23) = {p, j, m, A 2 , A 3 } (where 
p is the magnitude of the relative three-momentum of the 23 system, j and m are the 



C. Spectator equations for three spin 1/2 particles 



r 1 ) = 2M 1 G 1 v 12 \r 1 ) , 



(1.1) 



Y\) = 2Ml 2 G\V l2 \T\), 



(1.2) 
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fc x = -q 

(a) (b) 

FIG. 4. Diagrams showing the momenta in the two-body rest frame (left panel) and in the 
three-body rest frame (right panel). We chose the momenta of particle 1 to be in the — z direction, 
so the boost from the two to three-body rest frames is +z direction. 

angular momentum of the pair and its projection in the direction of q, and A2 and A3 are 
the helicities of particles 2 and 3, all defined in the rest frame of the 23 pair). The momenta 
are defined in Fig. 4, which also shows the relation between the rest frames of the two and 
three-body systems. The three particles have mass m, and the total mass of the three-body 
bound state is denoted by M t . [We use the symbol m to denote both the projection of the 
momentum and the particle mass, but the difference between them should be clear from 
the context.] Using this notation, and suppressing isospin, the final form of the three-body 
spectator equation for T 1 is given in Eq. (|6.3|). It can be written 

9crit 



(J1(23)p\T 1 ) = Y: E / ftflT d X sm X 

,//,//„,, Jo iy Jo 



x(j(23)p\Mi\j(2"3") P ")^g p "(q,f" 

tip" 

xPf7[l(2"3"), l'(2'3')] (J'l'(2'3')p|r 1 ) , (1.3) 

where (j(23)p\M 1 \j(2"3")p") is the jth two-body partial wave amplitude for the scattering 
of particles 2 and 3 in their own rest frame (precisely the amplitude obtained from the two 
body spectator theory as described in Ref. []2S]]), Vi 2 p [l(2 // 3 // ), l'(2'3')] is the matrix element 
of the permutation operator, given in Eq. (|1.6| ) below, and g p (q,p) the propagator of the 
off-shell particle in different p-spin states 

9+{9 >® = 2E^W 9 > g '^ P) = -W q ' (L4) 
where W q is the mass of the 23 pair, and depends on q, 

Wg = Ml + m 2 — 2M t E q , (1.5) 



with E q = ^Jm 2 + q 2 . Note that Eq. ( p..3|) includes a sum over intermediate helicities 
and angular momentum quantum numbers, and an integration over the internal spectator 
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momentum q' and the angle x between the directions of q' and q. The momenta p' and p" 
depend on q, q', and \i as given in Eq. (|5.17|) . 



The integration over q' has been limited to the finite interval [0, Qcrit], where q cr i t is the root 
of the equation W q ri = 0. At this critical spectator momentum (equal to ~ 4m/3 ~ 1200 
MeV), the two-body subsystem is recoiling at the speed of light and the relativistic effects 
are enormous! One consequence of this is that the solutions of the three-body equations 
go smoothly to zero as q — > q CT n (this is discussed in detail in Sec. VI). Contributions from 
q' > <?crit ; which come from two-body states with spacelike four- momenta, are suppressed 
both because of this zero and because the propagators for large q are small. Hence, even if 
the spacelike two-body scattering amplitude is not small, we expect spacelike contributions 
to the overall three-body amplitudes to be very much suppressed, and it seems sensible to 
simply neglect the region q' > g crit and set the three-body amplitudes to zero there. This 
also removes the need to calculate two-body amplitudes for spacelike total four-momenta. 

Exchanging particles 1 and 2 implies that particle 2 becomes the spectator and now 
its momentum and helicity must be expressed in the cm. frame of the three-body system, 
while the variables of particles 1 and 3 must be expressed in the rest frame of the 13 pair. 
Boosting from one frame to another introduces Wigner rotations of both the single particle 
and two-body helicities. In the helicity basis, this exchange operator is 



7< p '[l(2"3"), l / (2'3')] = (-I) m - Xl+X 3yj2j + Iyj2f + 1 

x 4iaj (#0 4^ (-#0 Nx»( 3 (q, q', x) , (1.6) 

where the functions d^ 2 ^ ((3) are the Wigner rotation matrices, and My, x , (q, q', x) describes 
exactly the Wigner rotations of the off-shell particle 3, as well as the nontrivial matrix 
elements between the different p-spinors u and v of particle 3 as they appear in the rest 
frames of the 23 pair and the 13 pair. The matrix j\f is defined in Eq. ( |6.2| ), the angles 9' 
and 9" in Eq. ( |5.17| ), and the Wigner rotation angles (3\ and /?2 in Eqs. ( |5.24| ) and ( |5.27[ ). 



For practical calculations it is more convenient to express Eq. (|1.3|) in terms of states with 
definite isospin and parity. These states will be denoted \Tj r {m\)p) , where we suppress 
reference to the total angular momentum and parity J n = l/2 + , T = or 1 and r = ±1 are 
the isospin and parity of the 23 pair, and A = A2 — A3 = 1/2 — A3. As discussed in Sec. VID, 
the states of good parity are superpositions of positive and negative helicity states [see 
Eq. (|6.17 )1, so that the two-body subspace is fully described by adopting the convention 



A2 = +1/2, and identifying the states by their parity r and helicity difference A = or 1. In 
this basis Eq. (|1.2j ) becomes 



(Tj r (mX)p\T 1 T )=Y: E / 9°^—/ d X sm X 



j'r' \" p" U ^1 
m'T' X' p' 

x (Tj r (m\)p\M 1T \Tj r (m\")p") -£L g»"(q,p") 

tip" 

xK P '[Tj r (m\")p", rfim'XW (T'j' r ' (m'X')^) , (1.7) 

tip' 
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where the permutation operator V P 12 P is given in Eqs. ( |6.27|) and (|6.30|) . Note that Eq. ( |1.7|) 
includes a sum over the intermediate isospin T' . 

This concludes our brief introduction; we now turn to a detailed derivation of the three- 
body equations ( p..3| ) and ( |1.7| ) given above. 

II. THREE-BODY EQUATIONS IN OPERATOR FORM 

We start with a derivation of Faddeev-type Bethe-Salpeter equations and introduce the 
spectator equations afterwards by substituting a new propagator and repeating the deriva- 
tion with all necessary modifications. 



A. Bethe-Salpeter Equations 

The total scattering amplitude for the three-nucleon system T can be decomposed into 
three parts T l , 

r = _r*. (2.i) 

i=i 

The partial amplitude T* sums up all diagrams in which particle i is the spectator during 
the "last" interaction (in the sense of "leftmost" in the diagrams of Fig. 1). Each amplitude 
T l is further split into sub-amplitudes T y , this time according to which particle does not 
participate in the "first" (or "rightmost") two-body interaction, 

3 

r* = J2 r ' lj ■ ( 2 - 2 ) 

The amplitudes T y satisfy the integral equation 

T ij = i 5 l] M l Gi 1 - M l G l BS T kj , (2-3) 

where Gi is the propagator of a single off-shell particle i, G l BS = G BS ®li = —i Gj^G^®li is 
the free two-body propagator for the {j, k} pair, and M. % = JkA 1 ®^ is the two-body scattering 
operator acting in the two-body subspace of particles j and k, with lj the identity operator 
for the spectator particle i. In our notation Gi is real, and any overall factor of i which 
emerges when the operator expressions are represented by Feynman diagrams is included in 
the propagator G BS . If V 1 = V l <S>h represents the sum of all irreducible diagrams describing 
the interaction of the two particles j and k with particle % a spectator, the Bethe-Salpeter 
equation 

M i = V { - V i G i BS M i (2.4) 
yields the scattering operator M. % . 
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A bound state of the three-body system can be defined as the residue of a pole of the 
three-body scattering amplitude T. For the triton we denote the position of the pole as 
P 2 = M 2 , where P = k% + k 2 + k 3 is the total four-momentum of the system and the ki are 
single-particle four-momenta. One can write T %: > as the sum of a pole term and a part 7£ u 



regular at P 2 = M 2 : 



Tl3 = -itf^ P L 2+' JZ ^ (2-5) 

where |P) are the partial vertex amplitudes for the bound state. Insertion into Eq. (| 
multiplication by (M 2 — P 2 ), and performing the limit P 2 — > M t 2 yields 

1") -VIY/^X F>. (2.6) 

These are the Bethe-Salpeter equations for the partial bound state vertex amplitudes. 

Up to this point the equations are very general and apply to systems of any three distin- 
guishable particles. Now we want to specialize to the case of three identical particles. We 
define the transpositions Vij of two particles % and j as follows 

Vi2 1 abc) = \bac) 

V 13 \abc) = \cba) , (2.7) 

Note that Vij interchanges the quantum numbers of the particles in the ith and jth locations 
in the state ket. The symmetry of the scattering amplitude under particle interchange can 
be expressed as 

V l3 T= CT 

TV l3 =(T } (2.8) 

where ( = +1 for bosons and —1 for fermions, and T is the symmetrized version of T. If 
we introduce the combined amplitude 

|r> = Ein, (2.9) 

i=l 

then the symmetry ( |2.8| ) of T carries over to i.e., 

^■|r>=C|r> 

(m,= C(r|. (2.io) 

These relations can be used to derive the permutation properties of the individual ver- 
tex factors If the particles are identical, then the two-body scattering operators and 
propagators acting in each two-body subspace are identical, and this is expressed formally 
by the relations 

Vij M i = M j 
VijG l BS Vij = G 3 BS , (2.11) 
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where M is a symmetrized version of M.. Using these, and the fact that V, 



1, we obtain 



—Vij M i Vij Vij G BS £ \T k ) 



Comparing with 



one obtains immediately 



C|r>> 



-M 3 G J BS Vij (\T) - |P) 
-M j G j BS (UV V,, V) 

= -m>g bs (cm-cn 



v lJ \r) = cm. 

Thus the three-body equations for identical particles can be written 

\V) = -(M i G BS (V lJ+ V lk )\r). 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



The three equations for the three possible choices of i are equivalent. It is therefore sufficient 
to solve Eq. ( |2.15| ) for, say, i — 1, and calculate |T 2 ) and |T 3 ) by means of Eq. ( |2.14| ). 

Eq. Q2.15D can be simplified further if we take into account the fact that the two-body 
amplitude M 1 is symmetric or antisymmetric under exchange of particles 2 and 3 for the 
case of identical bosons or fermions, respectively. Thus 



V23M 1 = M 1 V 2 3 = (M 1 . 

Using this relation, Eq. ( |2.15D (with i — 1) can be written 

T^ir 1 ) = -CPmMWbs (V 12 + Vrs) IT 1 } 
= -C 2 M 1 G 1 BS (P 12 + p 13 )|r 1 ) 

= Cir 1 ). 

Next, using the definitions Eq. ( |2.7| ) note that 

V23P12V22, 1 a6c}= V23V 12 \acb) = V 2 3\cab) 
= \cba) = Vi3\abc) . 

Hence, the operator "P 13 can be written 

*Pl3 — V23 r P\2V23 ■ 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



Using the relations ( |2.16| ) - ( |2.19| ) together with the fact that G BS commutes with V23 we 
can write the Faddeev equations (|2.15| ) in the following simple form 



ir 1 ) = -CM^ls (v 12 + ^23^12^23) ir 1 ) 

= -(M 1 G BS (l + (V 2 3)V 1 2\r l ) 

= -2CM 1 G^ s P 12 |r 1 ). 



(2.20) 



To reduce these equations to a practical form, it is sufficient to evaluate the permutation 
operator P 12 . 
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B. Spectator Equations 



Now we turn to the spectator equations. We begin by replacing the two-body propagator 
G BS ® lj, which describes the propagation of particles j and k (both not equal to i) in 
Eq. ( |2.6|) , by a new propagator, 

G i BS ®l i ^G j Q k ®l i , (2.21) 

where Q k is a projection operator which places particle k on the positive energy mass-shell, 
and, as in the BS case, Gj is the propagator of a single off-shell particle j. Choosing 
particle k to be the spectator during the "previous" interaction gives the unclosed form of 
the spectator Faddeev equations 

|P) = - ]T APGj-Qfclr*), (2.22) 

where the sum is over k and j with i fixed and no two indices equal. Explicitly, Eq. fl2.22p 
is shorthand for the following three equations: 

ir 1 ) = - {m 1 g 2 q 3 \t 3 } + M^sQsir 2 )} 
|r 2 ) = - {ii^GsQiir 1 } + M 2 G!Q 3 |r 3 )} 

|r 3 ) = - {M^iQair 2 ) + M 3 G 2 Q l \T 1 }} . (2.23) 

Note that the projection operator Q k insures that particle k is on-shell both as it leaves the 
partial amplitude \T k ) and as it enters the two-body scattering amplitude M\ 

To make a closed set of equations from Eq. ( |2.22D , it is necessary to place the final 
spectator particle i on shell, which then also forces one of the two interacting particles in 
the final state (denoted by k') to be on shell. The spectator scattering equations are shown 
diagrammatically in Fig. 2. The final bound state equations can be written algebraically in 
the following form 

QiQk'in = - J2 Qk>M l Q k Gj Q k Q t \T k ) , (2.24) 

where no summation over the index i is implied, and we used the projection property Q k Q k — 
Q k and GjQ k = Q k G r 

Alternatively, we may introduce the notation 

|r}> = D 

G{ = GjQ k <g> lj , (2.25) 

where the indices i, j, and k are all different, and the lower indices on M and T label which 
particles, apart from the spectator, are on mass shell. In this notation Eq. (|2.24j ) becomes 

\n>) = -Y, M l'kGl\T k ), (2.26) 
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We will use this notation in most of the remainder of this section, but will return to the 
definitions (|2T25|) later in the paper. As an example, consider the case i = 1 and k' = 2: 

\Tl) = -Ml 2 G\\T\) - M l 2Z G\\T\) . (2.27) 



As discussed in Ref. |J, f° r distinguishable particles without three-body forces Eq. (|2.26|) 



becomes a coupled set of six equations for the six amplitudes instead of only three 
equations for three as in the Bethe-Salpeter case. 

We emphasize that an important difference between the spectator subamplitudes 
and the Bethe-Salpeter subamplitudes |P) is that it is no longer possible to add the spectator 
subamplitudes together in order to construct a total amplitude, as we did in Eq. Q2.9|) . This 
is because the amplitude \T 2 ), for example, restricts particles 1 and 2 to the mass shell, 
while the amplitude \Tl) restricts particles 1 and 3 to the mass shell, and hence they are 
defined for different regions of phase space. Only operators or amplitudes which satisfy 
identical constraints, such as \Tl) and for example, can be combined. No total three- 
body amplitude exists in the spectator formalism. 

For identical particles, Eqs. ( ]2.26| ) can be further reduced by using permutation operators. 
Using the fact that the operator M 1 is symmetric under particle interchange, Eq. ( p,16|) , 
and the relation 

V 32 Q 2 = Q 3 V 32 , (2.28) 

we obtain 

V 32 M 2 \= CM\ 2 = Ml 3 V 32 

P 23 M 3 1 3 = CA4 = M\ 2 V 23 , (2.29) 

where ( = +1 for bosons and —1 for fermions, as before. (These are the operator form 
of the symmetry relations discussed in Ref. | 28fl ; note that, in the spectator formalism, the 



exchange operator does not relate an amplitude to itself, but to another amplitude with a 
different particle off-shell.) We will find it convenient to exploit the fact that Vjk = Vkj, and 
always write relations like those above so that the initial and final indices on both sides of the 
equation match. The spectator and on shell interacting particle can also be interchanged, 
leading to the following relations for the operators 

V\ 2 G X 2 V 2 \ = G\ 

V 23 GlP 32 = G\ . (2.30) 
The two-body amplitudes exhibit a similar symmetry 

V 12 M l 22 V 21 = M\ x 

V 23 M\^P 32 = M 22 . (2.31) 



Further relations can be found by combining the relations ( |2.29| ) and ( [2.31 ). One relation 
we will use below is 

V 12 M^P 32 V 21 = (Mf t = M* 3 V 31 . (2.32) 
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It is now easy to derive the effect of permutations on the spectator subamplitudes. For 
example, under the interchange of two particles in the interacting pair, 

V, 2 \Y\) = -V, 2 M l 22 G\\Y\) -V, 2 Ml z G\\T\) 
= -CMl 2 G\\Y\) - CMl,G\\Y\) 

= C\T\). (2.33) 
Using this, the interchange of the spectator with the on-shell particle in the interacting pair 



is 



V l2 \T\) = -V 12 M^V 21 V 12 GlV 21 V 21 \Vl)-Vi 2 M^Vs2r2iV 12 V 23 Gl\Vl) 
= -M^GlP 21 \Tl) - M 2 13 P 31 P l2 G\P 23 \Y\) 
= -M^GlP^Yl) - M 2 13 G 2 3 P 3l P l2 P 23 \Y\) 
= -M 2 u GjP 21 \Yl) - M 2 13 G 2 3 P 2l \Y\) 

= -M^GlP 21 \rl) - (M 2 3 G 2 3 \Y 3 2 ) , (2.34) 

where Eq. (|2.19|) (with 1 <->• 3) was used in the next to last step. Comparison with the 
equation for \Tf), 

C|r?> = -CM^Gl\Vl) - CM 2 l3 G 2 3 \Yl) , (2.35) 

implies 

P 2l \T\) = t\T\) . (2.36) 

Using these relations, we can obtain a single equation for IT^) 

\Y\) = -CMl 2 G\P 2l \Y\)-eMl 3 G\P 13 P 32 \Y\) 

= -CMl 2 G l 2 P 21 \Y\) - t 2 M l 23 P 32 P 23 G l 3 P 32 P 2l \T\) 

= -2CM 2 1 2 G^P 21 |ri>, (2.37) 



where Eq. ( 2.19 ) was used in the next to last step. From now on we will only consider 
fermions, so that the three-body equation for the vertex function T\, which singles out 
particle 1 as spectator in the "last" interaction and particle 2 as the interacting particle to 
be put on mass shell, becomes 

\Y\) = 2M 1 22 G\P 21 \Y\). (2.38) 

This equation is illustrated diagrammatically in Fig. 3. 

A more explicit form of Eq. Q2.38| ), which expresses M 22 and G\ as operators, and \Y\) 
clS db vector in Dirac space, is 

^2)^7 = 2[M 2 1 2 ] /3/3li77l [G , 2]/3 1 /3 2 , 7 i72 ^12 1^) q , (2.39) 

where a, (3, and 7 are Dirac indices for particles 1, 2, and 3 respectively, and summation 
over repeated Dirac indices is implied. Note that M and G operate on a two-body space 
only; the third particle (the spectator) is unaffected by these operators. 

In the next section we will give a momentum space representation of these equations. 
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III. MOMENTUM SPACE REPRESENTATION 



We specialize to three identical particles with mass m, spin 1/2, and four-momenta ki, 
k 2 , and k 3 . The total momentum 

P = k 1 + k 2 + k 3 (3.1) 

is conserved. The Gross equation restricts two of the three particles to be on mass shell, 
which for the choice Q2.38 ) are particles 1 and 2, with particle 3 off mass shell. In the 
three-body cm. system 

P = (M t ,0), (3.2) 

where M t is the mass of the three-body system, the momenta are 

ki = (£ fcl ,ki) 
k 2 = (E k , 2 ,k 2 ) 

k 3 = (k 30 , k 3 ) = (M t - E kl - E k2 , -ki - k 2 ) . (3.3) 

In Eq. (|3.3|) the four- momentum of the off-shell particle, k 3 , is fixed by four-momentum 
conservation. It is obvious that the problem has only 6 independent momentum variables, 
just as in the nonrelativistic case. 

The three-body basis states are direct products of a single particle state and a two- 
particle state, 

\h(k 2 k 3 )) = \h) ® \k 2 k 3 ) , (3.4) 

where, by convention, the off-shell particle has a bar over its momentum. 
Completeness and orthogonality relations are: 

(kx(k 2 k 3 )\k[{k' 2 V 3 ))= 2E kl S 3 (h - k[)2E k2 5 3 (k 2 - k' 2 )5\P - P') (3.5) 
1= / ^^d'PMhk^ihihk^l . (3.6) 

Next, we specify the matrix elements of all operators in this momentum space basis. The 
propagator is 

(k 1 (hh)\[Gl]i,0., r1 '\k[(k!M)) = 2E t ,S 3 (h - k[) 2E k2 S\k,_ - k' 2 ) S'(P - P') 

x2m[A + fe)] w , ^|^ , (3.7) 

where A±(k) = (m ± ^)/2m are the positive and negative energy projection operators. The 
two-body M matrix is 

(hfohMM^frwlWkgg)) = 2E kl 5\k 1 - k[) 8\P - P') 

xM^ ilY (k 2 3,k 23 ;P-k 1 ), (3.8) 
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where P — fa is the total two-body four-momentum, and the relative momenta in the two- 
body space are denoted by 

Note that kij = —kji. The two-body amplitudes in Eq. (|3.8| ) are identical to those discussed 



in Sec. IIA of Ref. [28]. The partial vertex amplitudes will be written 

{ki(fafa)\T\) al3l = r a py{kx, fa, fa) , (3.10) 

where, by convention, it is understood that the last momentum is the one which is off shell. 
Therefore 

{fa(k 2 fa)\Vi 2 \Yl) aPl = -(fa(fafa)\Tl)^ = -Tp a ^{fa, fa, fa) . (3.11) 
We can now obtain the momentum space representation of Eq. ( |2.38 ). Inserting the 



completeness relation ( |3.6| ) gives 

r r\ 3 V d 3 k' rl 3 k" rt 3 k" 

(fa{fafa)\T\) aPl = 2 / 1 ' ^(t.WI^I^CT 

k 1 «, 2 ""l /t 2 

MK(VM)\[Gl]ww^KWffi)) (k'((k^kl)\V 12 \Tl) a ^n . (3.12) 
Inserting the above expressions for M and G, and carrying out all integrals gives 



r a /3 7 (fci, fa, fa) — —2 d k 2 ——Mf 3/3 /^y( y k 23 , k 23 ; P — fa) 

J fcjU 

x [A+(^)]^ ( ? + ^ )Y7 : W(%, A$) , (3.13) 

772 /Cg 26: 

where fa = P — k' 2 — fa. This equation is manifestly covariant. 

These equations may be further reduced by multiplying the M matrix and the three- 
body vertex functions T by the on-shell spinors u (for on-shell particles in the initial state) 
and u (for on-shell particles in the final state): 

rA 1 A 27 (^i, fa, fa)= w a (ki, Ai)%(k 2 , X 2 )T a/3l (fa, fa, fa) 
M A 2 A^, 77 '( fc 23, fa 3 ; P - fa)= %(k 2 , X 2 )M 0P - m/ (k 23 , k' 23 ; P - fa)u^i\i! 2 , A 2 ) , (3.14) 

where w a (ki, Ai) is an on-shell Dirac spinor with three- momentum ki and helicity Ai. This 
gives us quantities with "mixed indices"; a Dirac index on a matrix element is replaced by 
a helicity index when it is contracted with a w-spinor of that helicity and with matching 
momentum. These amplitudes are still covariant, and simpler because the four- dimensional 
Dirac space is replaced by a two-dimensional helicity space. If we then replace the on-shell 
projection operator by a sum over on-shell w-spinors 

[A+(A&) V = E u ( k 2, A 2 )n(k 2 , A 2 ) , (3.15) 

and multiply Eq. ( |3.13j ) from the left by u a (kx, Ai)wg(k 2 , A 2 ) we get 
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/Tfl 
d s k' 2 —M X2X > 2m ,(k23, k' 23 ; P - ki) 



Equation ( |3.16| ) is still manifestly covariant, but is not suitable for a numerical solution. 
The main reason is that the two-body M matrix is given as a partial wave expansion in the 
two-body rest frame, and not the three-body cm. system, as needed in the above equation. 
A related problem is that the propagator for particle 3 depends on the angle between the 
vectors ki and k' 2 and is therefore not diagonal with respect to all angular momenta after a 
partial wave decomposition. 

In the nonrelativistic case, the first problem does not occur because the partial wave 
expansion is invariant under a Galilean boost, and the second is solved by introducing 
Jacobi coordinates. Because of the different energy-momentum relations in special relativity, 
neither of these problems can be handled so simply here. 

However, we can eliminate these problems here by exploiting the covariance of the for- 
malism, and by explicitly boosting the two-body subsystem to its rest frame. To prepare 
the way, introduce the total four-momentum of the two-body subsystem, 
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k 2 + k 3 = P -h = P + q, (3.17) 



where here it is convenient to introduce the momentum q = —k±. Next, the boost operator 
Afej is defined by the requirement 

A kl P 2 3 = P23 = (W q ,0). (3.18) 

The square of the mass of the (23) pair is then 

^ 2 = P 2 2 3 = P 2 2 3 = (M 4 -^) 2 -q 2 , 

= M 2 + m 2 - 2M t E q . (3. 19) 

A tilde ("~") on top of a variable always indicates that it is defined in the two-body rest 
frame. We have, e.g., 

^2 = A kl k 2 , 

h = A kl h- (3.20) 
We now define the relative momentum p through 

h = "^23 +p=(Ep,p) 

h = 1^23 -p=(W q -E p , -p) , (3.21) 



and therefore 



p=~k 23 = \{k 2 - k 3 ) = (Ep - \W q , p) . (3.22) 



Next, we introduce the representation on the Dirac space, S(A), of a Lorentz boost A. 
These transform Dirac matrices and spinors according to the following rules 
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S-\A)^S(A) = A» vl v 

S(A)u(k, A) =^P^ /2) (it! Afe ) W (Ak,/i) 



(k,A)5- 1 (A) = E^(Ak,/i)Pi 1 / 2) *(^ 



i-A/cJ 



(3.23) 
(3.24) 

(3.25) 



where -Ra& is the Wigner rotation accompanying the boost that connects the momenta k 
and Ak (see Appendix C), and the Dirac indices have been suppressed. The propagator of 
the off-shell particle 3 in the three-body rest frame can therefore be expressed in terms of 
its form in the two-body rest frame using ( |3.23| ) 



[m + 



m 2 — kn — ie 



s-'(A t .) <'"t fe) . s(A t ,; 



^3.26) 



^3 — m 2 — k 2 — ie 

Similarly, the full two-body M matrix in the three-body system can be written 

M 2>3 (k 23 ,k' 23 ;P 23 ) = S 2 \A kl )S^(A kl )M 2)3 (p,p';P 23 )S 2 (A kl )S 3 (A kl ), (3.27) 

where the subscripts 2 and 3 are shorthand for pairs of Dirac indices on particle 2 ((3f3', 
etc.) and on particle 3 (77', etc.), and the two-body scattering amplitude M(p,p'; P 23 ) is a 
solution of the two-body Gross equations in the two-body cm. frame. [Do not confuse M 2t3 
with amplitudes like M 23 used in the last subsection; here the subscripts refer to the Dirac 
indices, and in the previous subsection they referred to which of the interacting particles 
was on shell. From now on we have made the choice that particle 2 is on shell, and in the 
language of the previous subsection, all two-body amplitudes are M 22 .] Using ( |3.24j ) and 
( |3.25| ), we obtain the following expression for the mixed index M matrix 



M X2 v 3 (k 23 ,k' 23 ;P 23 ) 



S 3 \A k jV^(R Akik jM^ 3 (p,^ 



(3.28) 



where summation over all repeated indices (including helicities) is implied. Substituting 
these relations into Eq. (|3.16| ) gives 



F\ 1 \ 21 (ki, k 2 , k 3 ) 



-2 / d%^- ^(VJP^Ci^^M^^Cp^;^)^^^) 

k' 2 

(m + ^3)7273 



x- 



m 2 — k' 2 — it 



s 737 /(A fel )rv AlY (A4 h, k 3 ) . 



(3.29) 



This equation can be further reduced if we decompose of the propagator of the off-shell 
particle 3 into positive and negative energy parts 



(m + k 3i 



IT 



ni 



2 — k 2 — ie 



E 



E 



p a 3 
m x ^ 

~E~~ ^ 

a 3 



M 7 (k 3 , A 3 )M 7 /(k 3 , A 3 



-k 3 , A3)w 7 /(-k 3 , A 3 ) 



2E p -W q -ie 



-p, A 3 )« 7 '(-P, A3) V y (p, \ 3 )vy(f>, A3 



2E* - W n 



)€ 



W q -ie 



(3.30) 
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where the second expression can be obtained from the first using the fact that the spinors 
depend only on the three-momentum and k 3 = —p. At this point it is convenient to 
introduce "p"-spin by letting 

^■m:(-p a) a) .";=-. ^ 

Then the decomposition ( 3. 30]) becomes 

<(-P,A 3 )2 p (g,pK,(-p,A 3 ), (3.32) 



m 2 — k% — ie Ep 7 
where summation over p and A3 is implied, and 

1 



9 (q,p)-- 



2E P -W q -ie 



Substituting ( |3.32| ) into Eq. ( |3.25| ), multiplying from the left by u p {k 3 , A3), and using ( |3.24| ) 
and (|3.25| ), give the following reduced three-body equations 



x/(p',g)r^ A , A ,(^,^,^), (3.34) 

where R p Kk is the Wigner rotation for the spinor u p , and 

r AiA 2 A 3 (^i> k 2, h) = u p (k 3 , A 3 )r AlA27 (A; 1 , k 2 , k 3 ) 
KU^'P^'i W «) = u^^M^^p'; W^;04/4) . (3.35) 

We have reduced the three-body equations to six-dimensional integral equations for the 
coupled set of 2 4 = 16 amplitudes ^x 1 x 2 x 3 , which can be written 

Kx 2 x 3 (h, k 2 , h) = (k 1 X 1 (k 2 X 2 k 3 X 3 )p\r) . (3.36) 

The new states |fciAi(/c2A2/c3A 3 )/9) have simple completeness and orthogonality relations (de- 
veloped in the next section) which make them a useful starting point for further discussion. 

This form (|3.34j) for the three-body equations displays the Wigner rotations which appear 
when the two-body scattering amplitude is boosted from the overall three-body rest frame to 
its two-body rest frame. For practical calculations the equations will be further reduced by 
decomposing the amplitudes into partial waves, which will be discussed in the next section. 
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IV. ANGULAR MOMENTUM STATES 



In this section we follow the conventions of Wick M and define a basis of three-body 
partial-wave helicity states. Completeness and orthogonality relations are defined and the 
matrix elements of the propagator and two-body scattering amplitude are obtained. Using 
these states, the operator equations (|2.38|) are written directly in terms of the partial wave 
states. To obtain the final equations, the matrix elements of the permutation operator must 
be evaluated, and this is done in the following section. 



A. Construction of the states 

The three-body states are constructed in three stages. First, we construct the state of 
particle 2 and 3 in its rest system, choosing the momenta so that k 2 lies in the xz plane with 
k 2x positive, as shown in Fig. 4. By convention, particle three is off shell, and requires both 
u (p — +) and v (p — — ) spinors to describe its Dirac structure. This degree of freedom 
is referred to as the "p-spin" of the off-shell particle. Next, we boost the (23) system to a 
frame with three-momentum q = — k x in the positive z direction, and take the direct product 
of this state with the state of particle one with its three-momentum k x in the negative z 
direction. Finally, we obtain the partial wave states by an angular average over the Euler 
angles {$,9,0}, as defined below. In shorthand, this three-body state is denoted 1 1(23)) , 
to remind us that particles 2 and 3 are the pair which was boosted from their rest system. 

Begin with the construction of the state for particle 2 with momentum |k-2| = p pointing 
in the positive z-direction, and with helicity A 2 . This state will be denoted by \(p, 0, 0), A2), 
where the second two arguments in the parentheses are the polar and azimuthal angles of the 
momentum. The state with momentum pointing in an arbitrary direction can be obtained 
by applying a rotation operator Raq,^ = e~ l(t>Jz e~ iejy e~ llJz through Euler angles 0, 9, and 
7. For vectors without internal structure, we need only two angles, and following Wick 
characterize the states by the polar angles 9 and 4>, and represent the rotations by R^qq, so 
that 

|(p,fl^),A 2 > = i2 Mo |(p,0,0),A 2 >. (4.1) 
Note that this differs by a phase from the convention adopted in Jacob and Wick G9| and 



used in Ref. |pq| , where the rotation was defined to be R^^^ instead of R^qq- The phase 
difference is 

R^,-M °> °)> A 2> = ^ X2 R<f>,e,o\(P, 0, 0), A 2 > . (4.2) 

As discussed in Wick || , the new phase convention turns out to have significant advantages 
for the treatment of the three-body system, and gives identical results if (ft = 0, where the 
two-body states were previously defined . 



The state for particle 3, which in the rest system of the pair has a momentum of the 
same magnitude but opposite in direction, is defined 

I (P, 0, 0) , A3 p) = R^Iq I (p, 7T, 7r) , A 3 p) , (4.3) 
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where p = ± is the p-spin of the state (for more details, see the discussion below), and the 
following phase convention is incorporated into the definition of the state |(p, n, it), A3 p): 

I (p, vr, vr), A 3 p) = e-^R^m 0, 0), A 3 p) , (4.4) 

where S3 is the spin of particle 3 (in our case S3 = 1/2). The phase factor e~ l7TS3 is precisely 
what is needed for the definition Eq. (|4.3j ) to agree with the phase convention of Jacob and 
Wick [ SJJ for "particle 2" , which was used previously in Ref . |2J| . To see this, recall that 



R n ^ = e -i*J. e -i*J v = e -inj y( ,inj z ( 45 ) 

and hence 

e- mS3 R^,o\(P, 0, 0), A 3 p) = e-^e— -V^Kp, 0, 0), A 3 p) 

= e-^ S3 - A3) J Ro,.,o|(p,0,0),A 3 p) 
= (-l) S3 - A3 ^,o|(p,0,0),A 3 p), (4.6) 

as used in Ref. f28j. (According to our phase conventions, the value of J z is independent of 
p; see Eq. (A9) of Ref. ||28|| .) The two-particle state (in its rest system) is now written 

\(p, 9, 0), A 2 A 3 p)= R^oKp, 0, 0), A 2 A 3 p) 

= ^,8,0 (Ifo °' °)> A 2) ® \{P, vr, vr), A 3 p) ) , (4.7) 

where we emphasize that the phase e~ tnS3 is included in the definition of \(jp, it, it), A 3 p), as 
given in Eq. ( |4.4j ). Two particle states of definite total angular momentum and total helicity 
can be projected from these general two-particle states by integrating over the polar and 
azimuthal angles 



\pjm, A 2 A 3 p) = Vj d<t> r^sin^^ 2 _A 3 (0,^O)i?^ o |(p,O,O), A 2 A 3 p) , (4.8) 

•J "0 



2j + l\ 1/2 



r 2tv ptv 

d<f> \ 
10 Jo 

where we use the abbreviation 

* = 1=1^7 ' < 4 - 9 » 

The next step is to boost the two-particle state in the direction of the positive z-axis such 
that its total three-momentum becomes q. The required boost operator will be denoted Z q 
(it is equal to A^ 1 of the last section), and the (23) pair can be treated like an elementary 
particle with momentum q in the positive ^-direction, with "spin" j and "helicity" m, and a 
mass W g given by Wjj = (P — ki) 2 . The boosted state of the pair is no longer an eigenstate 
of the single-particle helicities, because a boost which is not in the direction of a particle's 
momentum mixes helicities. The three body helicity state is then constructed by taking a 
direct product of the boosted (23) pair and the state of the single particle 1 with a momentum 
of magnitude q in the negative z-direction. For consistency, the same phase convention is 
used to define the state of particle 1 that was used before to define particle 3, i.e.. 

\(q, vr, vr), A x > = e'^ R^Kq, 0, 0), Ai) . (4.10) 
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This gives a three-body helicity state with total 3-momentum zero and the momentum of 
the pair in the positive ^-direction: 

\(q, 0, 0),pjm, X 1 (X 2 X 3 )p) = |(g, 0, 0), Ai) ® \pjm, A 2 A 3 p) . (4.11) 

Following the convention for rotation of states first introduced in Eq. (P~T), the state in 
which the momentum of the pair is in an arbitrary direction is obtained from ( 4.11|) by 
applying the rotation ^$,0,0 

\{q, 6, $),pjm, X 1 (X 2 X 3 )p) = i2*,e, |(?, 0, 0),pjm, Ai(A 2 A 3 )p) ■ (4.12) 

Finally, the three-body helicity states with fixed total angular momentum J and projection 
M are obtained from the states ( |4.12|) by the angular average 

\qJM,pjm,X 1 (X 2 X3)p) = VjTd^ FdQ sin BV$^_ Xl ($, 0, 0) \(q, 9, $),p,jm, Ai(A 2 A 3 )p) 

= rjjm / d® dO sin 6 / d<f> d9 sin PV?*_ Al ($, 6, 0) 
J0 Jo Jo Jo ' 

x^Sa 2 -a 3 (0 ; £ 0)i? $ , e ,o {I (?, 7T, tt), Ai>®^fl Mo |(p, 0, 0), A 2 A 3 p)} , 

(4.13) 



where r/j is obtained from Eq. (|4.9| ) by replacing j with J . Note that this expression contains 
the two-body partial wave states (|4.8|) , and if we denote the rotation -R^^o by Ru, and 



dU = l^d® / r/OsmO 



^ff;-A 1 ( < f ) ©,o) = ^ ) ;_ Ai (f/) 

we have 

\qJM,pjm,X 1 {X 2 X 3 )p) ^r]jj dUV^>_ Xl {U)Ru {\(q, tt, tt), Xx)®Z q \pjm, A 2 A 3 p)} . (4.14) 

Another useful form of Eq. (|4.13j ) is obtained by exploiting the fact that a rotation about 
the z-axis commutes with a boost in ^-direction, so that the operation of the rotations on 
the (23) pair can be written 

R^&flZqR^fifl — R^,e,oRo,o,(j>ZqR Q Q = R$ t e,<t>ZqRo,e,o ■ (4-15) 
On the other hand, the rotation of particle 1 can be written 

R<s>,e,o = R<s>,e,4>Ro,o,-4> —> R$,e,(f>e %(t>Xl > (4.16) 

where the last step is obtained by letting i? ,o,-</> operate on \ (q, tt, tt), Ai), and recalling that 
this state is an eigenstate of J z with projection —Ax- Finally, noting that 

2>£U(^ ©, o)v%m -xM I o) = e^^g;_ Ai ($, e, ^^(o, e, o) 

= e iXl *V^_ Xl ($, 6, 0) d> mM _ X3 {6) (4.17) 
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shows that the factors of e lXl< ^ cancel, and that Eq. ( [4.1 3D can be written 



\qJM,pjm,X 1 (X 2 X 3 )p)=r ]j r ]] J dS V^_ Xl (S) £ d9sm 9 d% M _ X3 (6) 

xRslktXMXjqXjp), (4.18) 

where 

\K\{K*M A s)P> = I (<?, 7T, tt), Ax) ® ^i^-J (p, 0, 0), A 2 A 3 p) (4.19) 

is the three-body state in its canonical configuration in the xz plane with special 4-momenta 
k°, k%, and fcg, as shown in Fig. 4, i?,s = R$ } e t tf> is the rotation which carries the three body 
system from its canonical configuration to the most general orientation described by Euler 
angles 0, and 0, and 

/•27T /*7T /*27T 

dS = / rf$ / de sin 6 



JO Jo 

(5) = V^-m ($,©,0)- (4-20) 

Eq. ( |4.19| ) shows that the canonical three-body configuration is constructed by starting 
from a two-body state in the two-body rest frame where the relative momentum of the two 
particles is restricted to the xz plane with polar angle 9, then boosting this state in the 
positive z direction, and finally adding the spectator (particle 1) with momentum along the 
negative z-axis. Since the most general rotation R$ is performed after the boost Z q , the 
Wigner rotations that accompany the boost are all rotations about the y-axis, which greatly 
simplifies the calculation. 

The results ( |4.14| ) and (|4.18 ) are equivalent, and either may be used to evaluate matrix 
elements. 



B. Representation of the states 

In the previous subsection we showed how the states 

|fciAi(fc 2 A 2 fc 3 A3)p) = Rs \KK(.k° 2 KM\)p) , (4-21) 



introduced abstractly in Eq. ( p.36| ), are to be constructed. These states can also be written 



as a direct product of the momentum space plane wave states introduced in Eq. ( $M and 
Dirac helicity spinors: 

\k 1 X 1 (k 2 X 2 k 3 X 3 )p) = e-^ +s ^ \h{k 2 k 3 )) 

®R S ([Rx,«,o u(q, Ai)] a Z q R oAo {u p (p, A 2 ) [R^^o u p (p, A 3 )] 7 }) , (4.22) 

where a, /3, and 7 are the Dirac indices of particles 1,2, and 3, respectively, and all rotations 
are displayed explicitly, so that all spinor states in Eq. Q4.22j ) are "particle 1" states (in the 
sense of Jacob and Wick). Explicitly 
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u(p, A) 



with 



/ cosh(?]/2) 
v 2Asinh(r?/2) 



X(A) u(p,A) 



/ -2Asinh(r]/2) \ 
\ cosh(r]/2) / 



X(A) 



(4.23) 



and 







tanhr/ 



P 







EL 



(4.24) 



(4.25) 



Since the helicity spinors depend only on the magnitude of the three momentum p, we have 
used the notation v(p, A) = v(— p, A), so that the correspondence given in Eq. ( [3.3 1|) now 
becomes 



u{p, A) 
u(p, A) 



if p = + 
if P = - • 



(4.26) 



All of these conventions are consistent with our previous work (check Eq. (A9) of Ref. |28 
with % = 1). 

Using the representation ( |4.22| ), and the orthogonality relations ( |3.5| ) we obtain general- 



ized orthogonality relations for the three-particle helicity states Q4.21 ) 

(l'(2'3')p'|l(23)p) = (A ; iA / 1 (^A' 2 ^A' 3 )p , |A: 1 A 1 (A: 2 A 2 feA 3 )p) 

= (k[(k' 2 ¥ 3 )Mk 2 k 3 )) [u(q',X' l )R~l R- s }R s R 7T ^ u(q,X 1 ) 

x u(fi, A 2 ) R Q } /fi Z-, 1 R s }R s Z q R Q ~ fi u(p, \ 2 ) 

x ''f (fi, A' 3 ) R^ R Q } I Z^R- S } R s Z q R oAo R^ u> P (P, A 
= 2£ fcl 5 3 (£; 1 - k[)2E k2 5 3 (k 2 - k' 2 )8\P - P') 5 x[Xl 5^ 

x 



(4.27) 



Note that the states are not orthogonal in p space. Using Eq. fl4.23| ) gives 

u p \p, A' 3 ) u'(p, A 3 )] = <^a 3 O pV (p, A 3 ) , (4.28) 
where, if p = + is the first column and p = — the second, the matrix representation of O is 

o p > p (p,\)-- 



I 1 — 2A sinhr] ' 
I — 2 A sinh?] —1 y 

(r 3 ) p / p - 2Asinh^(ri) pV , 



(4.29) 



where sinh^ = p/m. It is also useful to express the covariant product 2Ek 2 8 3 (k 2 — k' 2 ) in the 
rest frame of the two-body subsystem, where k 2 = p depends on the polar angles 9 and 0. 
Hence the generalized orthogonality relation will be written 
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(l , (2'3 / )p'|l(23)p) = (k[X' 1 (k' 2 X' 2 k' 3 X' 3 )p'\k 1 X 1 (k 2 X 2 k 3 X 3 )p) 
= 8\' iXl 5y 2 x 2 5 X > 3 \ 3 O p/p {p, A 3 ) 

x2E kl 6 3 (h - K)2E p 6 3 (p' -p)5\P' - P) . 



(4.30) 



Even though the states Q4.21 ) are not orthogonal, the matrix 0(p, A) has a simple prop- 
erty which will enable us to carry out the calculation much as it they were. From Eq. (|4.29 ) 
we obtain 



[0(p, \)0(p, X)) plp = Vp cosh 2 r, = 6 P , P S . 
Using this we can show that the completeness relation for the states can be written 



(4.31) 



1 = Qa'a Q/3'l3 <5 7 '7 



d 3 ki d 3 p Til 



2Eu, 2E~ E- 



'-^d'P J2 |l(23)p')^(p,A 3 )(l(23)p| 



d 3 k x d 3 p m 2 4 x 



2E kl 2Ep Ep 



A1A2A3 p'p 

hXi(k 2 X 2 k 3 X 3 )p') O p , p (p, A 3 ) {k 1 X 1 (k 2 X 2 k 3 X 3 )p\ , ( 4 -32) 



A1A2A3 p'p 



where Qi = Q a ' a is the positive energy projection operator for particle 1 (with Dirac 
indices a and a') introduced in Sec. IIB. Eq. ( |4.32|) tells us that the states span only the 
positive energy sectors of particles 1 and 2 (which is sufficient) but they span the entire 
four- dimensional Dirac space for the off-shell particle 3. 

We will only describe the emergence of the factor 5y 7 in the derivation of the completeness 
relation ( |4.32|) . To see how this factor emerges, evaluate the sum over A3, p, and p' explicitly 
using (|Q3D 

9 

]Jl / 

~Ev H u p y {p,X 3 )O p/p {p,X 3 ) u 7 (p,A 3 ) 
a 3 p'p 

ffi 2 I / 

^2 { J2p u j'(P> ^)u^(p, A 3 ) -sinh 77 ^ 2 ^u p Y (p, X 3 )u?(p, X 3 



m 



,A 3 p 

1 —r 3 sinh 77 
t 3 sinh f] 1 

2 



A3, p^p' 



+ 



m 2 sinh 7] 



J 77 



sinh 7] t 3 
—t 3 sinh fj 



77 



5y 7 -=2 cosh 2 fj = <5y 7 . 
E~ 



(4.33) 



Subsequent operations by the rotations and boosts leave this factor invariant. In the same 
way, the sums over Ai and A 2 give the projection operators Qi and Q 2 . 

We will now use these relations to work out the generalized orthogonality and complete- 
ness relations for the partial wave amplitudes ( |4.18|) . 

C. Generalized orthogonality and completeness relations 



Using the generalized orthogonality relations (|4.30|) , the definition ( |4.21| ), and the nota- 



tion 
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J'l"(2"3» = \ q "J f M',p«fm>; \'((KK)p) , (4-34) 



we obtain 



(J'l'(2'3')p'\Jl(23)p) = (q'J'M',p'j'm'; X^X'^p'lqJM^jm; X 1 (X 2 X 3 )p) 

x f d9' sin~9> d% x ,(9') r^sinedg A2 _ A3 (0)(l'(2'3')p / |l(23)p) 
= 5 A;Al ^a 2 <^a 3 VJ'Vj'VJVj J dS' T>^P ml _ Xi (S') J dS T>$^_ Xl (S) 
x rde'sm9'd^ M _ X3 (§') dhmO d% M _ X3 (9) 

J « 

x0 pV fe A 3 ) 2£ fcl 5 3 (A; 1 - ^) 2£^ 3 (p' - p) 5\P' - P) . (4.35) 
Writing the 5 3 functions in polar coordinates 

2E kl 5 3 {k 1 - k[) 2Ep5 3 {p' -p) = 2E q -^-^- 8{cosQ' - cosG)^' - $) 

x2e /^ ~® 6(cqs6' - cos 0)5(0' - 0) , (4.36) 
pi 

allows us to integrate easily over dS' and c?#', giving S' = S and 9' = 9. The remaining 
integrals over dS 1 and <i# can then be easily done using the orthogonality properties of the 
T> and d functions. We obtain 

(J / l / (2'3')pV 1 (23)p) = 5j / j5 M 'M5j'j5 m/rn 5x' 1 x 1 5x' 2 x 2 5x' 3 x s O p/p (p, A 3 ) 

x2e / ( VV - ^7' ) 5\P> - P) . (4.37) 
q 2 p 2 

Using Eq. ( f4.32| ) and the orthogonality of the T> and d functions, the completeness relation 
for the partial wave states can be derived 



1 = Qa'a Q/3'/3 S-f'-y = f^^^d'P £ | ,71(23)//) O p , p (p, A3) ( J1(23) P | 

A 1 A 2 A 3 p'p 

= I ^-^^d 4 PY: \qJM,pjm;X 1 (X 2 X 3 )p')O p/p (p,X 3 )( q JM,pjm;X 1 (X 2 X 3 )p\ 



Jq iJ-^p J^p J MjrI 

X 1 X 2 X 3 p'p 



(4.38) 

Note that this is consistent with ( [4.37 ). 



D. Reduction of the equations 

The partial-wave expanded three-body equations can now be obtained directly from the 
operator equation ( |2.38j ). Restoring the projection operators Q, this equation is 
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Qi^r 1 = 2Q 2 M 1 G 3 Q 2 [QtQ 2 ] V 12 [QiQ 2 ] QiQ^ 1 , (4.39) 

where we have made frequent use of the fact that Qi operates only in the space of particle 
i, and hence commutes with all operators which operate only on particles j ^ i, and the 
property QiQi = Qi- Replacing the two terms [Q1Q2] <8> 1 with the completeness relation 
( |4.38| ), and using the relation 

Q1Q2I Jl(23)p) = I Jl(23)p) , (4.40) 

we obtain 

(Jl(23) lrH = 2 V V [ q" 2 dq"p" 2 dp" m 2 r q' 2 dq' p' 2 dp' m 2 

\ \ )r\ I / j / , / r) 771 0771 rp2 / OP OP p2 

J'M'j'm' A' 1 'A' 2 'A»p 4 P3 J ltj l" 2h P" ^P" J 2i V ltj V' ^P' 

A 'l A 2 A 3 p 2Pl 

x(Jl(23)p|M 1 G3|Jl"(2"3")p4)Op 4 p 3 (p ,/ ,A^) 

x(Jl"(2"3> 3 |P 12 |J'l'(2'3')p 2 > P2P1 (p', A3) (J'l'C^Oprir 1 ) . (4.41) 

This result anticipates the overall conservation of the total four-momentum, P, and the 
conservation of the quantum numbers J, M, j, and m by the operator M 1 G 3 . 

The matrix elements of the operator M X G 3 are most easily evaluated using Eq. ( [4.141 ). 
Recalling that M l G 3 operates only in the two-body subspace (23), the matrix element is 

(Jl(23)p\M 1 G 3 \J'l'(2'3')p 2 )0 P2Pl (p',X 3 ) 

= vjvji jdu J dn , v$^(U)v$^_ K {u') 

X ((<?, 71", 71"), \l\Ru l Ru>\{q> 7T, 7T), Al) 

x (pjm, A 2 A 3 pl^^-^/^a^^lp'j'W, A 2 A 3 p 2 >0 P2Pl (p', A' 3 ) . (4.42) 
The orthogonality relation for single particle states is 

(^Tr,^)^!!^ 1 ^!^,^,^)^!) = ((q,e,^),x 1 \(q , ,e , ,^),\' 1 ) 

= 5 XlK 2E q 5^(q-q') 

= 5 XlX ,2E q S ( q ~ q > 5 (cos® - cose')5($ - $') . (4.43) 

Inserting this, integrating over dll' (which fixes U' = U), noting that the boost operator 
Z q and the rotation operator Ru commute with M 1 G 3 , and carrying out the U integration 
using the orthogonality of the T> functions gives 

(Jl(23)p\M 1 G 3 \J'l'(2'3')p 2 ) P2Pl (p', A 3 ) 

S(q - q') 

= SjJ'dMM'Sjj'dmm'dxj^x'^Eq — 

x (j(23)p\M 1 G 3 \j(2'3')p 2 ) P2Pl (p', A' 3 ) , (4.44) 
where we have used the following shortened notation 

\j(2'3')p 2 ) = \p'jm,X' 2 X' 3 p 2 ) (4.45) 
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and removed from the two-body matrix element of M 1 G% the factor of 8jj>5 mm '. Now, 
including the projection operators, the propagator G\ in Dirac space is 



m 2 — — ie 



2m Q 2 il </(p', A3) /(g,p') </(p', A 3 ') 



(4.46) 



where we inserted the decomposition ( |3.32|) with g p (q,p) given in Eq. ( p.33| ). The only 
dependence of the two-body helicity states \j(2'3')p 2 ) on p 2 and A' 3 comes from the factor 
u p2 (p', A 3 ), which leads to the following property (where there is no sum over p') 



u»'(p',\£uO'(p',\'£ \j(2'3')p 2 ) = |j(2'3>') u p, (p',\'^(p',\' 3 



(4.47) 



Using Eqs. ( [4.46|) and ( |4.47|) , and recalling Eq. ( |4.31|) , the two-body matrix element in (|4.44|) 
is reduced as follows 

^(j(23)p|M 1 G 3 |j(2 / 3 / )p 2 ) P2Pl (p>, A 3 ) 

P2 



m 

= 2m—— 

Ep' 


E 0'(23)p|M 1 [^'(p',A 3 ')/(g,p') 


u p '(p'X 


)] |i(2'3 / )p 2 )C P2P1 (p / ,A 3 ) 


m 

= 2m—— 

Ep' 


E (i(23)p|M 1 |i(2'3")p , )/(?,p') 




^(p',A 3 ); 


o P2P1 (p',K) 


m 

= 2m—— 

Ep' 


E (j(23)p|M 1 |j(2'3")p / ) 9 P '(q,p') 5 X » X > O p , p2 

PV2A" 


(p',X' 3 )0 P2Pl (p',X' 3 ) 


Ef,' 

= 2m-^ 

m 


(j(23)p|M 1 |j(2 / 3')pi)^ 1 (g,p')- 






(4.4 



Inserting this result into Eq. ([4.41|) gives finally 

. m r p" 2 dp" m f p' 2 dp' m 2 

1 q uc 

J'M'j'm' A"A''p 3 



2 3 ' 
A 'l A 2 A 3 P2P1 



x(j(23)p\M 1 \j(2"3")p 3 )g^(q,p") 

x(Jl(2"3")p3^i 2 |J , l , (2 / 30p 2 )O P2Pl (p / ,A 3 ) (J'l^^SOpiir 1 ) 



In Appendix A we show that 



0'(23)p|M 1 | J (2"3")p 3 ) = 7§fe M A 2 ^,A3A»fe/^23) 



(27r) 3 m 



(4.49) 



(4.50) 



where AaA „(p, p"; P23) are the two-body amplitudes previously determined from Ref. [28 

Eq. (2.88). 2 ' 

To obtain the final three-body equations, the matrix elements of the permutation oper- 
ator must be evaluated, which will be done in the following section. 
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(a) (b) (c) 

FIG. 5. (a) The momenta k'",k 2 ,kf in their canonical configuration, (b) The canonical con- 
figuration of the momenta k"°, /c 2 °, k'^° which result from the interchange of particles 1 and 2. (c) 
Figure showing how /c 2 °, /C3 are rotated into fei,fc2>&3 by the rotation Ry which equals Ry 
only when ki, k 2 , k% line up precisely with k'f, k 2 °, k'£ . 

V. THE PERMUTATION OPERATOR V 12 

In this section we will derive the matrix elements of the operator P 12 which interchanges 
the states of particles 1 and 2. Using the shorthand notation given in Eq. (|4.34 ), the action 
of the permutation operator is 

Pia| Jl(23)p) = I J2(13)p) , (5.1) 

where the relation of the momenta of the individual particles to the relative momenta q 
and p, and to the quantum numbers j and m, is unambiguously determined by the order 
in which the single-particle names are written. For example, the state |J2(13)p) is one in 
which the second particle is the spectator with momentum q, and particles 1 and 3 are the 
pair with angular momentum j and m, and particle 1 has CM momentum variables p and 
9. More precisely, from the result ( 4.18Q , the state |Jl(23)p) is obtained by averaging over 



rotations Rs> of a state with the canonical configuration k'°, k 2 °, k!^ shown in Fig. 5a, and 
the state | J2(13)p) is obtained by averaging over rotations Rs" of a state with the canonical 
configuration k"°, k 2 °, k'^°, shown in Fig. 5b. Note that the two configurations are related 
by interchanging particles 1 and 2, but the definition of the states requires that k'" and 
k 2 ° both be in the negative z-direction, and that k 2 ° and k"° both lie in the xz-plane with 
positive x component. 

Any three-body configuration with canonical orientations in the xz-plane can be com- 
pletely characterized by three variables. For the vectors in Fig. 5a, these variables can be 
chosen to be |k' x | = q f , |k' 2 | = p f , and the angle x between k' x and k' 2 (recall that |k 2 | is 
the vector k 2 in the cm. of the pair). For the configuration in Fig. 5b, the corresponding 
variables are |k 2 | = q", \k"\ = p", and the same angle x- If the total cm. energy of the 
three-body system is fixed, then there is a constraint between these three variables, leaving 
only two independent. If q' and p' are specified, the angle % can be determined, and this 
was the approach taken by one of us previously [0]. However, the final equations are more 
tractable if q' and x are specified, and p' is determined by the constraints, and this is the 
approach we will take below. 
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Examination of the two configurations shown in Figs. 5a and 5b shows that the rotation 
Ry = R n ,x,o ( n °t to be confused with the rotation Ry> discussed below which carries the 
configuration shown in Fig. 5b into 5c) will bring them into alignment, or 

R v k'l°= k[° 

R Y k%°= k'° . (5.2) 

Furthermore, since the final momenta k\, hi, and A; 3 can be obtained either by rotating the 
(fc'°)'s through R s >, or the (/c //0 )'s through R s » (because they are equal), we have the relation 

ki = R s ,k'° = R S nk'l° = R s »RvK° ( 5 - 3 ) 

which implies 

R s > = Rs"Rv ■ (5.4) 
This rotation Ry will eventually emerge from the derivation below. 



A. Initial reduction of the matrix element 

We now turn to the details of the evaluation of the matrix element. Using the form 
Q4.18T) for the three-body state, the matrix element of V\i can be written 

(J'l'(2'3')p'\V 12 \Jl(23)p) = (q'J'M',p'fm\\[(\' 2 \' 3 )p'\Vu\qJM,pjmAi(X 2 \3)p) 

= (q'J f M>,p'3>m>, X' 1 (X' 2 X , 3 )p'\q"JM,p"jm, X 2 (X 1 X 3 )p) 

= VJ'VJVfVj J dS> J dS" V^? m ,_ K (S')V^* m _ X2 (S") 

x jT S> sin & jT d9" sin 9" d^_ x , {&)d^ Xi _ X3 0") 

x (A*Ai (k' 2 °X' 2 k$X' 3 )p'\ R~ S }R SII \K°UK°KW°^)P) ■ (5.5) 

Hence the matrix element depends only on the rotation R^} S n. This rotation will be equal 
to Ry after the constraints imposed by the evaluation of the ( k'° \ k"° ) matrix element 
have been realized, but until then this rotation will be denoted Ry = R a>x i^. Hence 
Rs" = Rs'R^ Rs" = Rs'Rv'i and using the group properties of the rotation matrices 

A 

The invariance of the group integration insures that 

J dS" = JdV , (5.7) 
and the orthogonality relation for the T> functions, 

/ dS'V^) ml _ x , {S ! )V { ^l{S') = 8j,j5 M , M 5 m ,_ KtK % , (5.8) 

J 1 V 7 
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allows the reduction of ( |5.5|) to 

(fl'(2'3')p'\V 12 \Jl(23)p} = 2n6j,j8 M , M Vfri J J dV V^, 



m— A2 



(V) 



d6> sin 9> jf d9" sin 9" d^ x ,_ x , (fl'^-Aa (*") 
x(k'°\[ (k'°X' 2 k*X' 3 )p'\ Ry \K°UK°KW \)P) • (5-9) 



,(i) 



We now define the vectors 



Ryk'[°= 


h 


Ryi k 2 °= 




Ry k 3 °= 


h 



(5.10) 

(where ki is not lined up along the negative z axis and equal to k'° until Ry = Ry). 
This rotation of the vectors k" into the vectors ki is represented in Fig. 5c. Guided by the 
discussion leading up to Eq. fl5~2]) and the representation ( |4.22| ) for the three-body states, 
the matrix element involving Ry is a product of a plane wave momentum space matrix 
element and Dirac space matrix elements 

(K°K (k' 2 °\' 2 -kJ\> 3 ) P rRy \K°\ 2 {K°\j4°\,) P ) 

= {k[\k' 2 k'i)\Ry\k^{k'{ k'r))xU 
= 2E kl 5^\k[° - k x )2E k2 5^{k' 2 ° - k 2 )5^{P' -P)U, 



(5.11) 



where 



U 



e^u(q', X[) R-^ RyZ q R oAo u(p, X, 



e - l7rS2 u(p\ X' 2 ) R-}, Z^RyR^ fi u(q, X 2 ) 
u p (p', A3) R~^ R~\, Q Z~, 1 R v , Z q R Q Q fi R 1T ^ ) QU fj '(p, A 3 



u$X uZ>M P = K(KK)P% A 2 (A 1 A 3 )p)) 



(5.12) 



will be referred to as the reduced matrix element of the permutation operator. Note that 
q = |kx| and q' = \k.[\ (as above). The matrix element (|5.12|) contains Wigner rotations 
which result from the fact that the helicities {A-} and {Aj} are defined in different frames. 
We first turn to the evaluation of the 5 functions on the right hand side of Eq. ( |5.11|) . 



B. Evaluation of the 5 functions 

In Appendix B it is shown that the two delta functions can be written 

6(p-Po)8 iP' ~ Po) 



2E kl 6^(k[° - k^Ek^ik'" - k 2 ) = 4E Po E plo 6(P)6(a - n 

x 5 (cos 6 — cos (9q) 5 (cos 6' — cos 9{ 



Po 



P'o 2 



(5.13) 
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with 



Po = po(q,q',x) p'o = Po(q',q,x) 
% = o Q (q,q',x) 0'o = do(q',q,x): 



(5.14) 



and 



\ 



Po(q,q',x) 

{Oo(q,q',x)} 



cos 



{M t - E q )Eg, + qq' cos x 
W q 

_ W q E q , - (M t - E q ) E Po(q 



q',x) 



qpo(q,q',x) 



The first two 5 functions insure that the rotation Ry = R, 



is now R 



7r,x',0) 



(5.15) 



and the delta 



functions in p and p' fix the angle x' to x- The angle x wm remain a variable, since we 
prefer to express the "allowed" magnitudes of the momenta p and p' as functions of x rather 
than the other way round. 

We now combine the expressions (|5.9|) , ( |5.11|) , and (|5.13|) and insert the result into the 
three-body equation ( |4.49| ). In doing this we must be careful to change the arguments 
of the matrix element (|5.11| ), which is expressed in terms of (p', q', A^A^Aglp, q, A1A2A3), to 
(p", q, AiAgA^lp', q', A^AgA^), so as to agree with the labeling used in Eq. ( [4.491) . Carrying 
out the dp' and dp integrations then gives 

(Jl(23) P |r 1 )= £ ^JTTv^TTT £ / q' 2 dq'^- [dx^X d%l Xl , m/ ^(x) 



2 3 1 
A 'l A 2 A 3 P2P1 



xO-(23)p|M 1 |i(2"3")p 3 } ^(^^e^^^-AgC^^-Ai^) 

9 

x 7^7^ <(P"' M^MpsI^ A^(AiA' 3 )p 2 )) Cpap! (p, A3) (J'l'(2'3')p 1 |r 1 ) (5.16) 



where 



p' = Po(q',q,x) p" =Po(q,q',x) 

0' = W,q,x) e" = 0a(q,q',x) 



(5.17) 



and ((p", Ai(A2A3)p3|p', A2(A' 1 A 3 )p2)) is the reduced matrix element defined in Eq. ( p.!2| ). 
This matrix element is calculated in the next subsection. 



C. Wigner rotations and the reduced matrix element 

It will be sufficient to define a Wigner rotation only for the special case when a spinor 
with helicity A and three-momentum in the right half of the xz plane is boosted in the 
positive ^-direction, as shown in Fig. 6. The boost is denoted by Z q [defined in Eq. QB4Q 1, 
the initial three-momentum of the state by p (with magnitude p and polar angle 9), the 
final three- momentum by p (with magnitude p = q' and polar angle 9 = ti — x), so that the 
Wigner rotation TZ(q, q', x) is defined by the relation 
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FIG. 6. (a) The canonical configuration of momentum for the calculation of the Wigner rota- 
tion. The helicity is +^ in this example, (b) The transformation of momentum and spin after the 
boost in the +z direction. The spin is now no longer aligned with the momentum, but rotated by 
angle with respect to it. 



Z q u(p, A) = Z q R Q g Lpu(0, A) = R 0j e,oL p TZ{q, q', x)«(0, A) , 



(5.18) 



where the representation of the pure boosts Lk (for k = p or p) in four- dimensional space- 
time is 



with 



/ cosh rj k 
V sinh i] k 

tanh rjk 



sinh rj k \ 



cosh rj k J 



(5.19) 



k 

Ek 



In Appendix C we show that 1Z is a pure rotation about the y-axis, 

fc(q, q', x) = Ro,p,o > 



(5.20) 



(5.21) 



and find the general equation for cos (3 as a function of q, q', and \- Since > (3 > n, [3 is 
uniquely determined by its cosine. Using the result ( |5.21| ), we have 



Z,u(p, A) = 5>(P, ") €l 2 \P) 



(5.22) 



We now are ready to evaluate each of the matrix elements in Eq. (|5.12|) , but we will make 
the substitution (p', q 1 , X^X^X'^p, q, A1A2A3) — > (p", q, XiX'^X'^p', q 1 , X'^^X^) so as to agree with 
the labeling used in Eq. ( |5.16|) . Noting that k\ = k'° implies that p' = q and 9' + x — 17 ( see 
Fig. 5) the matrix element for particle 1 becomes 



U 



(i) 



AiA' 



e™ si u 



e tnSl E A Roi,oRo, x ,oRo,8',ou(p', 

v 

e i ^J2[u(q,X 1 )u(q,v)\ ^(A) 

V 



(5.23) 
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where, using the function ft defined in Appendix C, Eq. ( |C14| ), 

fti = ft(q',q,X), (5-24) 
Similarly, to evaluate the matrix element for particle 2 use 9" = n — x, p" = q', and 



which gives 



l/ 2) 

A 2 A 2 



where 



Ra,ir,0 Ro,ir,—a 
Rn,a,0 Ro,—a,ir i 



e- l ^u(p", Ag) R-l lfi Z- 1 R v R K ^u{q', A' 2 ) 
e _ i7rS2 R-^ qR^ ^R^ ^ q u (q', X' 2 

V 
V 

e-^- 2 ^d^(ft 2 ) = e-^-^d^l(-ft 2 ) , 
ft2 = ft{q,q',X)- 



(5.25) 



d^(ft2) 



(5.26) 



(5.27) 



Calculation of the matrix element for particle 3 is complicated by the fact that its physical 
four-momentum is off-shell, while the four-momenta used in the definition of the u p spinors 
are on-shell. However, as shown in Eq. ( |3.31| ), the four- momentum of the negative energy 
spinor is identical to the four-momentum of the on-shell particle in the interacting pair, and 
an efficient way to proceed is to first express both of the u p spinors in terms of u~ . Then it 
will turn out that the boosts of both p-spin states can be evaluated in terms of quantities 
related to the on-shell particle in the interacting pair. 

To this end, note that u + can be expanded in terms of u~ and 7 5 u~ 



u + (p,X) = 7 5 



^ + 2A^ 7 5 
m m 



u (p, A) 



(5.28) 



where the matrix 



1 

1 



commutes with all rotations and boosts. Then, using the fact that 



R 



0,7T,0 



-i7T7 5 a2/2 



(5.29) 



where a 2 is the Dirac matrix, and using the explicit form of the spinors given in Eq. fl4.23Q , 
we obtain 



(5.30) 



Combining this with Eq. ( |5.28| ) gives a simple formula for R nj7T ,o u + (p, A). These two relations 
will be summarized 
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R,,n,ou p (p, A) = N P (X) u(p, -A) 
where u + = u is implied and 



N+(p,\) 



^+2A^7 5 



iV_(p,A)= 7 5 



m m 

These relations can now be used to evaluate the matrix element 



(5.31) 



(5.32) 



7/(3) 

u \'lp>,\' a p 



uP (P", A3) ^Trio^o }",o Z i lR -«,x,o Z q' R o,e>fl R -K,n,o u p {p', A3 
[u(f, -A 3 ') N p ,(f, X^R^^R^Z^R^^p', A 3 ) u(jT, -AC 
£ [ntf, -1/) i^ x>0 iV(p", AS) iW, A 3 )itV, Xl0 iV-xW? 

vv' 



(5.33) 



where iVp = 7 iVp 7 and, because we were able to write the states R^^p uP in terms of the 
positive energy on-shell spinors u using Eq. ( |5.31|) , the matrix elements of particle 3 have 
been expressed in terms of the Wigner rotations which already appeared in the treatment 
of particles 1 and 2. Since the N factors commute with the rotations, the matrix element 
can be further simplified as follows 



U \»P>,\> 3 P 



E ^ U 'd [ SUm u(q', -v') Ro,- xfl NAf, A 3 ') N p (p', A' 3 ) u(q, -v) 



7(1/2) 



= -£e^'^ A ^ [A p/p d 1 + B plp d 5 ] , 

vv' 

where the matrix elements d\ and d§ are 

[u(q', -V) i?o_ X) ow(g, -1/)] = cC^L^-x) rf i = ^f-i/(-x) ( c ' c ~ ^v'ps's) 
u(q', -u') iVx,o7 5 u(q, -v)] = ^_ V {- X ) 4 = &]_ V {- X ) {2i/s'c - 2uc's) 
with 



(5.34) 



c = cosh(7y g /2) 
s = sinh(ry g /2) 



c' = cosh (77^/2) 
s' = sinh (77^/2) , 



and 



sinh rj q 



m 



sinh 7] q i 



rn 



(5.35) 



(5.36) 



(5.37) 



From the explicit form of the N's, we obtain 



/ E p „Ep, s x/ , x , //'// 



Api p 



Bp'p 



\ 

( oV E p"P 



-4A 3 'A' 3 — 
, 6 m 2 

-2X'/- 



3 m 2 
m 



-2xf ^ 

-1 

-Ep" \ 
m 

/ 



(5.38) 
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The sum over v and v' in Eq. (|5.34|) can now be carried out if care is taken to remove 
all phases which depend on v or v' . There are four possibilities, all of which occur. We may 
write the "standard" sum in a compact form 



E^ = E 



W ,(1/2) fa w (l/2) f _^Jl/2) 
e a -u>-\>. 



■(ft)<t'^_„(-x)^'4;(A) 



vv' 

JttX'AI/2), 



K ' 3d xl\l(fo + fa - X) = e™Mfa +P2-X) 



(5.39) 



where symmetry properties of the d functions have been used, and for compactness the A3A3 
indices are suppressed in the final result. By similar arguments, the remaining three sums 
give 



]T4z/VJ = 

vv' 




E2^ = 

vv' 


= e inS3 


E 2 ^ = 

vv' 


= e inS3 



d{-i3 l + (3 2 - X ). 



Using these identities, the matrix element for particle 3 finally becomes 



U 



(3) 



yP3p2 



w 3 ' 



where the matrix X p p is 
/ 



VP P 

A 3 A 3 



[DxA + D 5 B] 



'_ 1 U/2-Aj 



^_l)l/2-Ai 

with the notation 



V 



E ff pf 
m m 



D i = 4[S 3 (A + & - X)c'c - ^ (A + fa + X )s's 
A = A 



—As 

m 



m 



/(V2), 



++ 

B = 2A 3 5 +H 



(5.40) 



(5.41) 



(5.42) 



(5.43) 



Combining the results (|5.23[) , (|5.26|) , and (|5.41|) gives the following expression for the 
reduced matrix element 



AiA' x KPV d \>>\> 2 \~PV <*\»\' s 



P3P2 



(5.44) 
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D. Symmetries of the permutation operator 

The reduced matrix element (|5.44j) satisfies a symmetry condition which can be obtained 
from the following property of scalar products 

(if, X 1 (X'^p 3 \p', X' 2 (KX^p 2 )) = ((P\ m'^piW, A 1 (AM/> 3 )}* ■ (5.45) 

Because the permutation operator is hermitian and the initial and final states are composed 
of identical nucleons, this equation tells us that the matrix elements ( |5.9| ) must be identical 
under the substitution q <-> q' (which also implies p" <-> p', 9" <-> 9', and (3\ <-> /3z) and 

j / Ai A' 2 

m <r+ m! X 2 <r-¥ . (5.46) 

P3 P2 A^' A3 

Examination of the matrix elements shows that this implies 



4)-- 2 -3^ A , m _ Ai ( x )^(^)^(- A )^, (5.47) 
which reduces to the condition 

(-l)V^'^ = ^| . (5.48) 

However, the transformation (|5.46|) gives A <->• A, <-» — (— 1) A 3 _A 3_B, £>! <-» (— l) A s'— ^3Z? l5 
and D 5 <-> — D 5 , showing that the result (|5.42|) satisfies the symmetry condition ( |5.48| ). 

Another symmetry of the matrix elements of the permutation operator follows from the 
fact that V12 commutes with the parity operator V, which leads to the identity 

Vn = V 12 V 2 = VV 12 V . (5.49) 

The action of the parity operator on the states defined in Eq. Ql.l3| ) can be worked out, 
giving 

V\qJM,pjm,X 1 (X 2 X 3 )p) = {-iy-\-l) J - 3 - ll2 P \qJM,Pj -m, -A x (-A 2 -A 3 )p) . (5.50) 

Hence the matrix elements of Vi 2 must satisfy the identity 

(q' JM',p'j'm', X[(X' 2 X' 3 )p'\V 12 \qJM,p 3 m, X 1 (X 2 X 3 )p) 

= (q' JM',p'j'm', X' 1 {X' 2 X' 3 )p'\VV 12 V\q JM,pjm, X 1 (X 2 X 3 )p) 

= (-l) 2J - 1 pp'(q'JM',p'f -m^-X^-X'z-X'Jp'lVulqJM^j -m, -Ai(-A 2 -A 3 )p) . (5.51) 

For the triton, where J = 1/2, this means that under the substitutions 

X[ <-> —X[ X\ <-> —Ax 

X' 2 <— > — X' 2 X 2 <— > — A 2 

A3 <-> — A3 A3 <-> —A3 

m' <-> — m' m <-> — m (5.52) 
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we should recover the same matrix element multiplied by a factor of pp' . 

To verify that the matrix elements of Vyi satisfy this symmetry, return to the full ex- 
pression given in Eq. (|1.6|), and use the identity 

d%i§) = (-l) X - X 'd U {_Ad) (5-53) 

to obtain the condition 

A/jg; = (-1) A 3-VVA/-S ' , (5.54) 



where Afy/y is defined in Eq. (|6.2| ) below. Examination of this equation confirms that 
Eq. ( |5.54|) is indeed satisfied. 

In the next section we present our final results for the three body equations. 



VI. FINAL EQUATIONS 

In this final section we collect the previous results together, and explain how it is that 
integration over the the spectator momentum, q', is limited to a finite interval. Then we 
describe the changes in the equations which are required by isospin and the conservation of 
parity. Finally, we describe how the three-body channels are classified and counted. 



A. Spectator equations in angular momentum space 

Using Eqs. (gjp and (fO|) gives the following compact result 

Yl x \l>tpp*pi(p '> ^3) : 



P2 



E-i 

m A 3 A a 



where the matrix J\f p ' p is 



\fpp 

A 3 A 3 



in 



3/V3 ^5 

m 



V 



-2X' 3 D 5 



-2A 3 [DiB + D 5 A] 



Eff p' 



m 



m 



J 



and A, B, Di and D 5 were defined in Eq. G5.43| ). 

Combining Eqs. ( |5.44j ) and ( |6.1| ) and substituting into Eq. ( |5.16| ) gives 



(ji(23)pir i )=^ v /2j+i v / 2 J '+i x: 



A~ A„ p 



9crxt m 

Q dq — 

E q i JO 



(6.1) 



(6.2) 



dXSinx^-Ai.m'-A'M 



'12 3 ' 



x( J (23)p|M 1 | J (2"3")p")^/'(g,/)rf 0) 



xf-1 



,m-Ai+A' j(1/2). 



Epii 



l")d%,„ x ,{0') 



d^(-P2)fV^^ (J'l'(2>3')p>n . (6.3) 
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This is identical to the final result given in Sec. I, Eqs. (|1.3|) and (|1.6|) . It is a two-dimensional 
integral equation depending on the variables q' and Xi where the integration over the angle \ 
runs from to tt, independent of the value of q', and the integration over q' has been limited 
to the finite interval [0,g crit ], as discussed in the next subsection. The momenta p' and p", 
the angles & and 9", and the Wigner rotation angles /?i and fli all depend on q, q', and Xi 
and are defined in Eqs. (|5.17|), (|5.24|) and (|5.27|). The matrix M is defined in Eq. 



B. Removal of the space-like region 

The physical reason for restricting the q' integration in Eq. ( |6.3D to the finite interval 
< q' < (fa-it will be discussed now. 

As given in Eq. Q3.19I ), the invariant mass of the two-body subsystem decreases with 
increasing momentum of the spectator, q, and at the value 

M? — m 2 4 . , 

Q — <7crit = — -Tn 6.4 

H H 2M t 3 V ; 

the mass of the two-body subsystem is zero. This means it is recoiling with the speed of light, 
and under such circumstances the relativistic effects are clearly enormous! Furthermore, as 
q increases beyond the critical value, we pass from a region where the two-body states 
are time-like into a region where they are space-like. The two-body scattering calculations 
are carried out in the rest frame of the two-body system, which does not exist for space- 
like states, and, more generally, it is unlikely that an effective theory designed to describe 
time-like scattering would be useful in the space-like region. Furthermore, since the space- 
like two-body states appear only at rather high momentum (above 1200 MeV) where the 
amplitudes are expected to be very small anyway, it would be sensible to simply neglect 
the region q > g crit , and set the three-body amplitudes to zero in this region. As it turns 
out, the three-body amplitudes go to zero automatically at the critical value ofq, permitting 
us to impose the condition that they are zero for q > q CI i t without making the three body 
amplitudes discontinuous in q. 

To see that the Faddeev amplitudes ( Jl(23)p|r 1 ) — > as q — > g crit , note that the function 
Po(q, Q' ' i X) [defined in Eq. ( |5.15| )1 approaches infinity as q — > q cr i t (as long as q' > 0, which is 
true over the entire region of the q' integration except at the boundary where the integrand 
is zero). Specifically, 

~n _ gcrit (E q , + q'cosx) C 

p ™ t w q ~wv [b - b) 



and in this limit 



m . 1 



g-(j>")N-»^ -± (^) K_ p -> {W q f , (6.6) 



E p ,, ^ ' W q \mCJ 



where C and K± p are functions which are finite in the limit as q — > g crit . Note that the 
possible 1/Wg singularity from the negative energy part of the propagator is canceled by the 
m/Epir factor. Hence the Faddeev amplitudes go like 
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(ji(23) P |r 1 ) Ti — - c p+ (wg^+ +1 ) + Cp _ (w q ) np - , (6-7) 



where n p+ and n p _ are powers with which the two-body amplitudes (f4.50|) fall with momen- 
tum as ft" — > oo: 

(1 \ n PP3 
^ • ( 6 - 8 ) 
V I 

We conclude that the Faddeev amplitudes not only go to zero as q — > g crit , but that they 
approach this limit smoothly. 

C. Isospin 

Since the main application of the three-body equations for spin 1/2 particles will be the 
three-nucleon system, we have to incorporate the isospin degree of freedom. This can be 
done separately from the other degrees of freedom, as described in this subsection. We will 
assume that isospin is conserved by the equations. 

To lay the foundation we return to the discussion in Sec. II. The exchange operators 
Vij are a product of a part which acts only in isospin space, and a part which acts on all 
other coordinates, denoted by V%j. If the ij pair has isospin Ty, the action of Vij on the 
isospin part of the wave function will be denoted simply by its eigenvalue (— l) TiJ_1 . The 
phase ( which occurs in Eq. ( |2.29| ) is —1 for fermions and is therefore a product of the phase 



(— l) TiJ 1 from the exchange of the isospin variables, and the phase u, resulting from the 
operation of Vij. Hence 

C = -l = n(-l)^- 1 . (6.9) 

Even though ( is always —1 for fermions, there are in general two possible values of u 
corresponding to the two possible isospin channels, and Eq. ( |2.29|) generalizes to 

V 23 M^= u M23 = M' 2 V 23 ■ (6. 10) 



The vectors \Tl) are also vectors in isospin space. Taking matrix elements of Eq. ( p. 38 ). 
and inserting 1 = J2t gives 

(T\Tl) =2J2(T\V 12 \T') M£G\Vi 2 {T'\Y\) , (6.11) 

where |T) are the isospin wave functions discussed below, and (T\V\ 2 \T') is the matrix 
element of the permutation operator in isospin space. The calculation of this matrix element 
is familiar from the nonrelativistic theory, but for completeness we will briefly present it here. 
In more detail, the states in isospin space will be denoted 

\T) = \{{t 2 h)Tt 1 )TT z ), (6.12) 
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where U is the isospin of particle i, T is the isospin of the pair, and T and T z are the total 
three-body isospin and its projection. As the notation suggests, t 2 and t 3 are first coupled 
to T, and then T and t\ are coupled to T. These states form a complete, orthonormal basis 

(((t 2 t3)Ttorr,|((44)T^)r%'} = 5 tlt ^ t2t ^ t 3 t ^ TT ^rr'^ 

J2 |((t 2 t3)Tt 1 )TT 2 )(((t 2 t 3 )Tt 1 )TT 2 | = 1. (6.13) 
The effect of V\ 2 is to interchange particles 1 and 2: 

Pi2\((t 2 t 3 )T , t 1 )TT g ) = |((tit 3 )r't 2 )rr,) . (6.14) 

The matrix element of V\ 2 in isospin space reduces therefore to a simple recoupling coeffi- 
cient, 

(T\V 12 \T f ) = (((t 2 t 3 )Tt 1 )TT 2 |P 12 |((t 2 t 3 )T , t 1 )Tr 2 ) 
= (((t 2 t 3 )Tt 1 )TT z |((t 1 t 3 )T / t 2 )TT 2 ) 

= -v / 2TTTv / 2T 7 Tl| ^ %t>}- (6 ' 15) 

In the next subsection we complete the reduction of Eq. ( |6.11| ) by inserting a complete 
set of good parity eigenstates. 



D. Parity eigenstates 

The two-body scattering amplitudes which drive the three-body equations are separated 
into channels which are eigenstates of V 1 (the parity operator on the 23 subspace) and 
isospin. Isospin was just discussed in the previous subsection, and need not be revisited 
again until the next subsection below where we explain how isospin (or exchange symmetry) 
plays a role in the description and counting of the channels. The role of the conservation of 
parity, which has not yet been taken into account, will be discussed in this subsection. 

The helicity states | Jl(23)p) are neither eigenstates of the full parity operator, V, nor of 
the two-body parity operator, V . Since the two-body scattering amplitudes which emerge 
from the two-body calculations are eigenstates of V , the three-body equations must be 
re-expressed in terms of these states. This is not difficult because the eigenstates are merely 
linear combinations of the states we have already obtained. 

First we return to the two-body helicity states |jm(A 2 A 3 )p), where the relative momen- 
tum, p is suppressed because it will play no role in the discussion which follows. If we apply 
the operator V 1 (referred to simply as V in Ref. P5| ) to this state we get 

P 1 |jm(A 2 A 3 )p) = pe \jm{-\ 2 -A 3 )p) , (6.16) 

where e = (— l)- 7-1 . It is easy to see that the state 

I j r (m\)p) = -L (l + rV 1 ) \jm(X 2 X 3 )p) 

= i= {|jm(A 2 A 3 )p> + rpe\jm{-\ 2 -A 3 )p)} , (6.17) 
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with A = A2 — A3, is a normalized eigenstate of V 1 

V l \j r {m\)p) = r \j r (m\)p} . (6.18) 

If we replace the individual particle helicities A2 and A3 on the rhs Eq. ( 6.17| ) by — A2 
and —A3, we obtain the same state, apart from a phase factor. We should therefore include 
in our new basis ( |6.17D only states that are not related to each other by changing the sign of 
both helicities. We choose the convention A 2 = +§ and label the states by the difference A. 
With this convention A can be or 1, and the parity r can be + and — , and we have again 
4 independent states, just as before when each of the individual helicities were allowed to 
be ±~. The selection rule 

|A 2 -A 3 |<j (6.19) 

excludes A = 1 for states with j = 0. 

For the construction of three-body parity eigenstates we can proceed in precisely the 
same way, treating the two-body subsystem as one elementary particle with spin j, helicity 
m, and (now well defined) intrinsic parity r. The parity operator which acts in the three- 
body space will be denoted V, and should be distinguished from V 1 . To carry out this 
construction, we first introduce the three-body states 

\Jj r \ 1 (m\)p) = -j= {\J,jm; Ai(A 2 A 3 )p) + rpe\J,jm; Ai(-A 2 -A 3 )p}} 

= J- (1 + rV A \JJm; X 1 (X 2 X 3 )p) , (6.20) 



V2 

where | J, jm; Ai(A 2 A 3 )p) = | Jl(23)p) are the same three body states introduced in Sec. IV, 
Eq. ( [4.34| ), but with some of the notation restored for clarity. The parity operation V on 
these states yields 

V\Jj r X 1 (mX)p) =r ll r{-l) J -^\J ] r -Ai(-mA)p) , (6.21) 
and the three-body eigenstates of parity, with eigenvalue IT = ± are therefore 

I J u j r (m\)p) = -L {\Jj r A!(mA)p) + n^r(-l) J -^|^ T -Ai(-mA)p)} 

= -L(l + TlV)\Jj r \ l (m\)p}. (6.22) 

In this case we adopt the convention Ai = +| and let m vary, subject to the condition that 

\m — \i\ — \m — -\ < J . (6.23) 

Nucleon one is always in a positive-energy state, and therefore 771 = 1 and s± = ~. The 
triton is characterized by J n = | + . Three-body states with these quantum numbers are 

j r (m\)p) = -L {\y \{m\)p) - re\y r -|(-mA)p)} . (6.24) 



V2 
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Expanding these three-body parity eigenstates in terms of the original three-body helicity 
states Q4.34D gives 

\\ + j r {m\)p) = - {||, jm; Ai(A 2 A 3 )p) + rpe|§, jm; Ai(-A 2 -A 3 )p) 



-re||,j -m; -Ai(A 2 A 3 )p) - p\\, j -m; -Ai(-A 2 -A 3 )p)} 
1(1 + p) (l + r^H^mjA^Aa)^. 



(6.25) 



We now return to Eq. ( |6.11|) and take the remaining spin-momentum matrix elements of 
the operators using the basis of good parity states Eq. fl6.25p . For simplicity, we represent 



the states by |~ j r (mA)p) = \j r (mA)p), and the direct product states by |T)(g)|j r (m\)p) 
\Tj r (m\)p) . Then, in our abbreviated notation, 



(Tj r (mX)p\T 1 T ) = E / ftf^- d X 



smx 



j'r' \" p" 
m'T' A'p' 



m 



x(Tj r (mX)p\M^ \Tj r (m\")p")—g<> (q,f) 

xV^ , [Tj r (my')p",T'f r \m'X')p'] ^ (T' j" {m'X')p'\T x T ) , (6.26) 

tip* 

where {Tj r (mX)p\M 1T \Tj r (mX")p") = M^b, (Tj r ) is the two-body scattering amplitude 
for the jth partial wave with parity r and isospin T, and 

v( 2 p '[Tj r (mX")p",T'j' r '(m'X')p') = (T\V 12 \T') x (j^(mX")p"\V 12 \f r ' '(m'X')p') 

= (T\V l2 \T') x ((Via)). (6.27) 

Eq. ( |6.26| ) is our final result. It expresses the three-body equations in terms of the physical 
states with definite parity and isospin, driven by two-body amplitudes M^b, (Tj r ) which 
have been previously calculated as described in Ref. . 

The new matrix element ((V12)) is readily obtained from the original matrix elements 
Eq. ( |5.5|) , which are, in the notation of this section, 

«Pi2» = (h j'm'; KMK)A'Pu\im Ai(A 2 A 3 )p) . (6.28) 
From the definition ( |6.25 ) we have 

((VT2)) = lihfm'^'MX^p'l (l+r'V 1 ) (1+V)V 12 (1+V) (l+rV 1 ) ||, jm; X 1 (X 2 X 3 )p) 



%j'm'; A^A^A^p'l (1 + r'V 1 ) V 12 {1 + V){1 + rV 1 ) ||, jm; X 1 (X 2 X 3 )p) , (6.29) 



where we used the fact that V commutes with V\ 2 . Using Eqs. (|6.16 ), ( 6.21 ), and ( |5.51|) , 
we obtain 
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((Vn)) = i{(i,i , m';A / 1 (A^)p'|P 12 ||,im;A 1 (A 2 A 3 )p) 

+rVe'(|, jV; A' x (-A 2 -A 3 )p'|7> 12 |§, jm; A 1 (A 2 A 3 )p) 
+rpe(l,j'm'; X'^X'^p'lV^Jm; Ai(-A 2 -A 3 )p) 
+rVp'pe'e(i,j"m / ;A / 1 (-A 2 -A 3 )p'|P 12 |i, jm; A x (-A 2 -A 3 )p) 
-pd^'W^UA^Ayp'i^ll^- -m;-Ax(-A 2 -A 3 )p) 
-r'p'pe'^j'm'; A;(-A' 2 -A 3 )p'|P 12 ||,j -m;-A x (-A 2 -A 3 )p) 
-re(~, jW; A^AaA^p'lP^H, j -m; -A x (A 2 A 3 )p) 

-r'rp'e'ed, jW; A' x (-A 2 -A 3 )p'|P 12 |±, j -m; -A 1 (A 2 A 3 )p)| , (6.30) 

This is the correct form of the permutation operator to be used with the physical, good 
parity states. 



E. Three-body channels 

We conclude this paper by counting and classifying the channels which contribute to the 
final three-body equations. 

In order to clarify the following discussion we restore some of the notation which we 
previously suppressed, and denote the two-body helicity states ( |6.17| ) with good parity, 
\j r (mX)p), by |p j r (mX)p), where 



p = ^-o-^23 <0.3U 



is the difference in the energies of the two particles in the two-body rest frame (with particle 
2 on shell) and jy 23 the rest frame energy of the two-body system. Note that po is in general 
not zero because particle three is off-shell, and that we continue to suppress explicit reference 
to the magnitude of the three component of the relative momentum, p, because it will play 
no role in the discussion which follows. This state satisfies the relation 

Q 2 |poj r (mA)p) = |p j r (mA)p), (6.32) 

where Q 2 is the projection operator introduced in Sec. II which places particle 2 on shell, 
and is equivalent to the identity when operating on states where particle 2 is already on 
shell. Note that the state with relative energy — po has particle 3 on shell and hence 

Q.3 I -Po j r (™A)p) = | -p j r {m\)p) . (6.33) 

In terms of the states (|6.17|) , the interchange of space and spin coordinates (everything 
but isospin) has the following effect: 

V 23 \PoJ r (mX)p) = (rp) x e^ | -p j r (m\)p) . (6.34) 

This equation follows from the definition (|6.17|) and the relation 
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^23 \Po jm{\ 2 Xz)p) = e | -Pojm(\ ?J \ 2 )p) , (6.35) 



which (except for notational changes) is Eq. (2.97) of Ref. |28| . 



Using Eq. ( |6.34|) , we can extend the discussion of the previous subsection and introduce 



states with both good parity and good exchange symmetry. Introduce the states 

\j™( mX ) p )==J-(l + uV 23 ) \p j r (m\)p), (6.36) 



y/2 

These are normalized eigenstates of both V 1 and V23 

V 1 \j ru (mX)p) = r\j ru (m\)p) 

V 23 \j ru (mX)p)=u\j ru (mX)p). (6.37) 

Since u = (— 1) T , which can be written T = (1 — u)/2, these states are also the correct 
spin-momentum states to use with isospin. 

The counting and classifying of three-body states depends in part on the number and 
classification of the two-body scattering states. Since both V 1 and V23 are conserved by 
the two-body equations, two-body states can be classified by different possible values of the 
quantum numbers r and u. There are four combinations: 

singlet r = —e u — e 

triplet r = — e u = — e 

coupled r = e u = e 

virtual r = e u = — e . 



The last set of states, referred to as virtual states in Ref. ||28|| , do not contribute to physical 
two-body scattering. This is because, in the positive p-spin sector, their parity assignment 
would require that j = £ ± 1, which in turn requires a total spin S — 1. These assignments 
are consistent with an exchange symmetry of u = — e only if these states are odd under 
change of sign of the relative energy variable po, which insures that they are zero on shell. 
However, because the two-body quantum numbers r and u are not conserved in three-body 
scattering, they can contribute to relativistic three-body scattering and to the three-body 
bound state. In the calculations completed thus far we have neglected these states, but 
we expect them to give a small contribution of purely relativistic origin. 

Neglecting the virtual states, and recalling the selection rule (|6.19|) leads to the following 
counting rules: 

j = : (A = 0) x (p = ±1) x (r = ±1) x (u = -1) = 4 states 

j > : (A = 0, 1) x (p = ±1) x (r = ±e) x (u = e) 

+ (A = 0, 1) x (p = ±1) x (r = -e) x (u = -e) = 12 states 

The total number of two-body states with angular momenta j < j max is therefore n 2 = 

4 ~\~ 12 jmax- 

In numerical calculations of the three-nucleon bound state it has become customary 
to truncate the partial wave series according to the maximal included total pair angular 
momentum j. Table | shows how many different three-body states exist for a given j. The 
pattern is simple: applying the selection rule ( |6.23j ) for each j > gives 24 possible states 
corresponding to 12 two-body states with either m = or 1, or 2 x 12 = 24 states. For j = 
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TABLE I. Possible quantum numbers for three-body states with J n = 5 (for the triton). 



Virtual two body states have been neglected. 



j 


A 


r 


u 


P 


rn 


number of states 








± 




± 





4 


> 1 


0,1 


± 


e 


db 


0,1 


16 


> 1 


0,1 


— e 


— e 


± 


0,1 


8 



only m = is allowed, and hence there are only 4 different states. For each combination 
of quantum numbers there is one particular pair isospin consistent with the Pauli principle. 
Since we have used exchange symmetry to count the states, the inclusion of isospin does not 
lead to any further increase in the number of channels. The total number of states up to a 
given maximal value j max is therefore n.3 = 4 + 24j max . 
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APPENDIX A: OPERATOR FORM OF THE TWO-BODY EQUATIONS 

In this Appendix we present the operator form of the two-body equations, and their sub- 
sequent reduction to partial waves. Using the notation of Sec. IIT3, the two-body equations 
for the scattering amplitude are 

M = V -VG2Q1M , (Al) 

where V is the symmetrized kernel and M is the two-body scattering amplitude describing 
the scattering of particles 1 and 2. Note that particle 1 in on-shell in the intermediate 
state, in agreement with the conventions of Ref. [f2SJ (see Sec. IIB in that reference). In the 



three-body language, the M in Eq. ( |A~l"| ) is M 3 , because the "spectator" (if it were present) 
would be particle 3. To obtain a closed set of equations for (|A1|), multiply by <2i, giving 

[QiM] = [Q X V\ - [Q1VG2Q1) [QiM] , (A2) 

which shows that particle 1 is on-shell throughout the interaction. 

The two-body partial wave equations can be obtained from Eq. (|A^) by inserting a 
complete set of the two-body angular momentum states defined in Eq. ( |4.8| ) [with the 23 
pair relabeled 12]. The completeness and generalized orthogonality relations for the two- 
body states, implied by the work in Sec. IIC, is 
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(j'(1'2')p'|j(12)p> = S fj 5 m >m5 K xA' 2 M O plp {p, A 2 ) 2E p S(p ^ 5\P' l2 - P 12 ) 

l = Q«> a 5p>p = j P ^%d±PY,\m)p')O p , p {M(m)p\, (A3) 

V P jm p 1 p 

AjA 2 A 3 

where the shorthand notation defined in Eq. ( |4.45|) has been used. Substituting the com- 
pleteness relation into Eq. (|A2"D gives 

(j(12)p|M|j(l'2')p'> = (j(l2)p\V\j(l'2')p') 

-E / ^^(jWVm^>i)0 P1P2 (KX^ (j(l"2> 2 |M|j(l'2')p') . (A4) 

P1P2 J "^fc 
x"x" 
1 2 

The evaluation of the matrix element of VG is identical to the evaluation of M l G% carried 
out in Eq. ( ggg ). Substituting this (|A|) gives 



(j(12)p\M\j(l'2')p') = {j(12)p\V\j(l'2')p') 



2 

m 



^ /fc 2 c^MJXl2)p|V|JXl' / 2 , Op , 0^/'(0,A;)(J(l' / 2' / )p' / |M| J (l'2 / )p , )• (A5) 

p" A'/Al,' -^fc 

Multiplying both sides of the equation by (2n) 3 m 2 / EpEpi , and introducing the amplitudes 
^^(j(12)p\M\j(l'2') P '}= M*l X2X , (p,p f ; P 12 ) 

(27r)3r? ^0-(i2)p|v|j(i / 2')p / >= ^aTa^a^^p;^) 

k ~ g p (0,k)=gp(k), (A6) 



EpE p 



(2tt) 3 
gives the equation 

" E / ^^A 2A „(p, fc; Pu)gf{k)MM^{k,fr Pn) ■ (A7) 
p" A"A^' 

This is identical to Eq. (2.88) of Ref. [2S|, establishing the relationship between this paper 
and previous work on the two-body problem. 

APPENDIX B: REDUCTION OF THE 5 FUNCTIONS IN V 12 

In this Appendix we derive Eq. (|5.13|) for the delta functions 

2E kl 5®(k' 1 ° - h)2E k2 5^(k' 2 ° - k 2 ) (Bl) 
which appear in the matrix element of the permutation operator, Eq. ( |5.11| ). 
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1. Radial 5 functions 



First find the vectors k'° and k' 2 ° shown in Fig. 5a, and k'[° and k 2 ° shown in Fig. 5b. 
The vectors k° and k 2 ° are 



(E q ,\ 



\-q'J 



k 2 ° 



( E i\ 




\-qJ 



:b2) 



where we anticipate that {k"} — > {k} and hence use [Ic^ I = Q instead of |k^'| = q" . This 
agrees with the notation in Eq. ( |5.12| ). The vectors k 2 ° and k x can be found by boosting 
the vectors k 2 ° and k x ° defined in the two-body rest frame: 



k 2 



( E p , \ 
p' sin 9' 


y p' cos 9' J 



/ E r \ 
p" sin 9" 


\ p" cos 9" ) 



(B3) 



The boost in the z-direction for the configuration shown in Fig. 5a is 



(C S'\ 

0100' 

10 

V s 1 c J 



(B4) 



where 



a 



Eq' 



Wl + q' 2 



S' 



The boost for the configuration shown in Fig. 5b is obtained from 
q. Hence 



(B5) 

by replacing q' by 



k 2 



Z q ,k° 



k\ — Z q k l 



(C'Ep + S'p^ cos 9' \ 
p' sin 9' 


\S'Ep, + C p' cos 9' ) 

( CEp,, + Sp" cos 9" \ 
p" sm9" 


\SE § ,, + Cp" cos 9" ) 



(B6) 



Now the vectors which appear in the delta functions ([Bl"|) are k'", k' 2 °, ki, and \s.2- The 
second two of these are obtained by applying the rotation Ry to k 1 ° and k 2 °, as illustrated 
in Fig. 5c. Since the rotation Ry does not change the length of the three-vectors, the 
condition k 2 ° = k 2 imposed by one of the radial delta functions becomes 
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q = {p' 2 sin 2 9' + S a El + C' 2 p' 2 cos 2 9' + 2C'S'p'E i> , cos 0'} 

1/2 



1/2 



C'E P , + S'p' cosd' 



— m 



Requiring that | cos#'| < 1 gives the unique solution 



cos 9' = cos 9' n 



E q - C'Ep _ W q ,E q - {M t - E q ,)E p 



S'p' 



q'p' 



The radial delta function can therefore be rewritten 



8(k' 2 °-k 2 ) =5(cos9' -cos9' ) 



dko 



The derivative is 

giving finally 

S(k' 2 ° - fc 2 ) 

A 9 



dko 



d cos 9' 



cos y'=cos t 



d cos 9' 

S'E q p' 
Q 



-i 



COS P'=COS t 



- \- 5(cos9' - cosO'n) = W * S(cos9' - cos^') , 
S'qEgp' V 07 qq'E q p' V 0,1 



(B7) 



vB8) 



(B9) 



where cos^q was defined in Eq. (|B~8|). 

The result for the other radial delta function follows immediately by interchanging primed 
and unprimed variables. 



2. Angular 5 functions 



Evaluation of the angular delta functions requires explicit consideration of the rotation 
Ry = R a ,x',/3 which rotates the vectors k x ° and k 2 ° into k\ and k 2 as discussed in Sec. VA 
and illustrated in Fig. 5c. First we consider the delta function 

5 {2) (k'° - k 2 ) = <5(cos 9 2 ° - cos 9 2 )5{<p' 2 ° - 2 ) , (BIO) 

where (9 2 ,(j) 2 ) and {9 2 ° , 4> 2 °) are the polar and azimuthal angles of k 2 and k 2 °, respectively. 
The three- vector k 2 ° was given in Eq. (p6|) , and k 2 is 





( 






( sin x' cos a \ 


k 2 — R a>x i t pk 2 ° — R a ,x',p 




: .) 


= -q 


sin x' sin a 










\ cos x' ) 



Since k 2y = k' 2y = and k 2x = k' 2x > 0, the azimuthal part of Eq. (B10|) becomes 

5(0 2 O - <fo) = S(n - a) . (B12) 
The polar part of the delta function is 
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5(cos02° - cos9 2 ) = 5 Q [S'Ep, + C'p'cosO'] + cosx'j • 



(B13) 



If this delta function is used to eliminate the x'i integration, one is left with an integration 
over p' with upper and lower limits depending on the two external momenta and on the 
second integration variable. This makes numerical solutions of the resulting equations awk- 
ward. Because x' is the angle between the momenta q and q', and is therefore symmetric 
under interchange of initial and final states, it is more convenient to retain x' as the indepen- 
dent variable and eliminate instead the integration over p' . The limits on the x' integration 
turn out to be independent of the other variables, running over the expected range from 
to n. Using Eq. ([BHD to replace cos 9', the delta function becomes 



5 (cos #2° — cos 82 



( C'E q — Ef 



+ cc&x 



with the solution 



Ep> = E plQ 



C E q + S'q cos x 



S'q 

(M t -E q ,)E q + qq' cos X 



(B14) 



(B15) 



and, since the k 2 ° = k 2 is now satisfied, we have replaced x' by X-> as discussed in Sec. VA. 
It is easy to show that Eq. ( B14|) implies that Ep> > m for all values of Xi so that the delta 
function ( |B14j ) places no additional restrictions on the x integration. Hence 



^(cos^ — cos 62) = S(p — p'q) 



d cos 6*2° 



dp' 



-1 



p'=p'o 



S'qE &o 



P'o 



S(P - P'o) ■ 



The final angular delta function is 

5 {2 \ti° - h) = 6{cos6[ - cos^ 1 )5(0' 1 ° - 00 



(B16) 



(B17) 



where (0i,(pi) and (9°,(f)°) are the polar and azimuthal angles of k x and k x °, respectively. 
The vector \i° is given in Eq. ( |B2| ); the vector ki is 



kl — -Rvr^^k]^ — Rtt,x,P 



( Vx\ 




' v x cos (3 \ 1 







v x s\n(3 = - 


w 




\ v z J \ 



v x cos (3 cos x + v z sin x 

v x sin f3 
v x cos (3 sin x — v z cos x 



(B18) 



where 



v x = p" sin 6" 



v z = SE r +Cp" cos 6'' 
Setting ki = k'" gives three equations 



CEqi — Epii 

S 



(B19) 
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0= v x cos (3 cos x + v z sin \ = sin #1 cos 4>i 

0= f x sin (3 = q' sin 6*i sin 0! 

g= v x cos /3 sin x — v z cos x = —q' cos #i . 



(B20) 



We will first use the left hand set of Eqs. ( |B20|) to obtain the values of (3 and p" which are 
fixed by the delta functions ( B17 ). Then we will use the right hand set to find the Jacobian 
of the transformation from the variables cos 82, 02 to the variables p" , (3. 

The angle (3 must be or 7r. Allowing for either possibility, the first and third of 
Eqs. (|E20D give 



q cos x = _ v z . 
Substituting this result back into the third equation gives 

q' sin 2 x — v x cos (3 sin x = P" sin 9" cos (3 sin x ■ 



(B21) 



(B22) 



Since the sin of all angles under consideration is positive, this equation shows that cos (3 is 
also and that therefore (3 = 0. Finally, from Eq. ( |B21| ) we obtain 



cosx + 



CE 1 — Ef; 



q'S 







(B23) 



which is the condition (P14j ), with q and q' interchanged, showing that Epn satisfies ( |B15|) 
(with q and q' interchanged). We have shown that the second angular delta function is 



SW(k' 1 °-k 1 )=6(p)6(f-p )J : 



(B24) 



where po is p' with q and q' interchanged and J is the jacobian of the transformation from 
the variables cos 6 2) 4>2 to the variables p", [3. 

We return to the right hand set of Eqs. ( |B2U| ) to calculate this jacobian. Unlike the 
previous cases it is necessary to calculate a full jacobian because the variables are all coupled 
unless we go to the limit sin^ = 0, which gives singular results. Postponing this limit until 
the end, we first eliminate v z from the Eqs. ( |B20| ) and obtain 



v x sin (3= q' sin 9\ sin $1 

v x cos (3= q' [sin 6\ cos 4>\ cos x ~ cos &i sin x] ■ 

Differentiating both of these equations with respect to cos 9\ and 
and then taking the limits Q\ = n, 0i = 0, and (3 = gives 



dv x 


J = 


dp" 





d<t>i 
dp 



d cos 81 d cos 8\ 



q 



cos x ■ 



(B25) 

>i, computing the jacobian, 
(B26) 



The radial delta functions have fixed cos 9" in terms of p" — > po [Eq. ( |B8| ) with q' — * q and 



p — > p 



po], and using this relation we find that 



dv x 
dp" 



q'W q cos x 
qEp sin 9" 



(B27) 
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and hence the jacobian is 



J 



gg'-Epo = Sq'E. 

PO W q po 



PO 



(B28) 



Combining Eq. flBlf), its companion, Eqs. ( |B16| ), ( |B24| ) and (B28), and anticipating the 
fact that the delta functions fix the rotation so that p" — > p and 6" — > 0, we obtain our final 
result 



Po 



x 6 (cos 9' — cos9' )5(cos9 — cos# ) • 



(B29) 



APPENDIX C: WIGNER ROTATIONS 



In this Appendix the Wigner rotation angle for the standard boost ( 5.18|) that occurs in 
the matrix elements of the permutation operator is derived. 

Consider a spin 1/2 particle with mass m, helicity A, and three-momentum p which lies 
in the xz plane (with the x-component positive, by convention). Under the boost Z q in the 
+z direction the state transforms like 

Z q \p, A) = R(Q)\p, A) = £ \p, MZJPiP) » ( C1 ) 



where Z q will by used to denote both the boost ( p4]) in four-dimensional space-time, and its 
representation on the space of states. Hence p = Z q p. In agreement with the notation used 
in Sec. IV, the magnitude |p| = q', and angle between p and the +z axis is 9 = n — \. We 
will show that Q is a rotation about the y axis, and find the rotation angle, (3, in terms of 
q, q', and X- 

As previously described, the helicity states are constructed from rest states by first 
boosting in the +z direction, and then rotating through the proper angle. For the states 
\p, A) and \p, A), this construction gives 

\p,\)=e- iJ » S L p \m,\) 

\p,\)=e- iJ » e L ql \m,\), (C2) 

where L k was defined in Eq. (|5.19| ) and rh = (m, 0) is the four-momentum vector of a particle 
of mass m at rest. Hence the Wigner rotation operator is given by 

R(Q) = (e- iJ y L q ) ^ Z q e~ iJ y' e L p . (C3) 

Assuming that R(Q) is a pure rotation about the y-axis, the equation reduces to the following 
set of equations involving (5 

e~ iJy6 L q , e- iJy(B = Z q e- iJy§ L p (C4) 
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We will solve these equations using the Dirac representation for the operators. In this 
representation the pure boosts in the ^-direction are 



c + s a z 



(C5) 



where tanhr/g/ = q'/E q >, d = cosh (j] q i/2), and s' = sinh(^ (? //2) [these relations were previ- 
ously defined in Eq. ([5736])]. The boost Z q has a structure similar to ( |C5| ) but with r] q / — > n q 
where 



C 



M t -E q 



cosh n n 



S 



sinh n n 



(C6) 



as in Eq. (|B5|), but with q' — ► q. To simplify notation in what follows, we denote the 
hyperbolic functions of n q /2 by 



c n = cosh(n 9 /2) 
The pure rotations about the y-axis are 



s n = sinh(ra 9 /2) 



e -u y e 



-ie^a y /2 = cos (#/ 2 ) - i sin(^/2) 7 5 a ?; . 



(C7) 



Hence Eq. ( P4|) becomes 

cos(#/2) - i7 5 Oy sin(0/2)] [c' + s'a 2 

= [c n + a z s r 

Using {7 5 a y ,a 2 } = 0, Eq. ( |C9|) becomes 

{cos[(6 + /3)/2] - rfoty sin[(0 + P)/2] } d + {cos[(# - /3)/2] - i^a y sin[(0 - p)/2] } s'a z 



cos(/3/2) -i7 5 a y sin(/5/2) 

cos(#/2) - i7 5 a„ sin(0/2) 



[cp + a 2 s p ] . (C9) 



\C n Cp -\- S n Sp + a z (s n c p + c n s p )} cos(6 l /2) 

-* {c n c p - s n s p + a z (s n c p - c n s p )} 7 5 a y sin(0/2) . 



(CIO) 



Equating the coefficients of the independent operators on each side of this equation gives 
four coupled equations 

d cos[(# + P)/2] = (c n c p + s n s p ) cos(^/2) 

d sm[(0 + P)/2] = (c n c p - s n s p ) sm(6/2) 

s' sin[(0 - P)/2] = (c n s p - s n Cp) sin(0/2) 

s' cos[(0 - p)/2] = (c n s p + s n c p ) cos(0/2) . (Cll) 



These are not all independent; squaring each of these and then adding the first two and 
subtracting the second two gives an identity (1=1). Hence only three are independent, and 
given the quantities q, q', and x — 71 — Q, these three independent equations can be regarded 
as equations for the unknown quantities 0, p, and (3. 

We are only interested in an equation for j3. This is obtained by multiplying the first 
equation by the last equation, and adding it to the product of the second equation and the 
third. The result is: 
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C smh Tjp -\- S cosh Tj p cos 9 

cosp = — . (^12) 

smh 7]' 

To eliminate cos 9, we find an equation for it by summing the squares of each of the equations. 
This gives 

~ cosh rj — C cosh ri v 

cos9 = — — p - 

b smn rjp 

= ^t. (C18) 

Note that this is identical to Eq. ( |BE| ) [with the primed and unprimed variables exchanged] , 
showing that the calculations are consistent. Substituting for cos 9 gives 

cosh r] p cosh -n' — C 

cos ^ = — r r — r ~r^ — 

smn T] p sinn rf 
W q E p E q ,-m 2 (M t -E q ) 



pq'W q 

q'(M t -E q ) + qE q , cos X 



lq' 2 W 2 + q 2 E 2 q , + 2qq'E q ,{M t - E q ) cos X + {qq' cos X Y 
This is the formula for cos \]3(q, q', x)]- 



(C14) 
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